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ABSTRACT 


From  experimental  pressure  drop  and  flas?  rate  data  of  the 
pilot  pipeline  tests,  the  apparent  viscosities  of  various 
crude  oils,  for  the  temperature  range  of  20  to  S0°F.,  were 
calculated.  Possible  viscosity  dependency  on  shes-r  rate, 
as  for  a  non-Newtonian  fluid,  was  investigated  by  vanriation 
in  flow  rates. 

The  relationship  between  the  calculated  apparent  viscosity 
and  temperature  was  defined  by  a  single  curve  for  a  specific 
crude  for  both  laminar  and  turbulent  flow  pattems  regardless 
of  the  apparent  shear  rate  at  the  wall,  the  diameter  of  the 
test  line  or  the  storage  tank  temperature  of  the  crude  oil. 
Apparent  viscosity  data  for  both  turbulent  and  laminar  flow 
conditions  were  found  to  be  identical  in  the  relatively  wax-free 
range  of  40  to  80°F.  The  crude  oils  showed  linear  apparent 
viscosity- temperature  relationship  for  oil  temperatures  in  the 
range  80  to  65°F.  Belov/  55°F  the  apparent  viscosity  of  all 
crudes  increased  rapidly  with  decrease  in  temperature.  The 
transition  occurred  in  the  range  65  to  55°F. 

Even  under  laminar  conditions  of  flow,  the  pilot  pipeline 
investigations  showed  the  flow  of  crude  oil  through  the  test 
lines  to  be  characterized  by  high  apparent  shear  rates  at  the 
wall.  The  shear  rates  observed  may  have  been  large  enough 
to  give  the  limiting  viscosity  characteristic  of  a  pseudoplastic 


% 


fluid  at  infinite  shear.  The  pseudoplastic  nature  of  similar 
crudes  at  temperatures  below  their  "wax  points"  where 
precipitated  wax  is  held  in  suspension  mthin  the  fluid, 
has  been  reported  by  previous  investigators.  The  maximum 
shear  rates  used  by  these  investigators  were  appro:djnately 
40  reciprocal  seconds  while  minimum  shear  rates  at  the  wall, 
obtained  in  the  pilot  pipeline  tests,  were  rarely  below  1000 
reciprocal  seconds  and  exceeded  100,000  reciprocal  seconds 
at  the  higher  rates  of  flo?/.  The  results  indicate  that, 
within  the  operating  range  of  the  pilot  pipeline,  the 
crudes  tested  behave  as  Newtonian  liquids.  However, 
pseudoplasticity  cannot  entirely  be  disregarded  in  view 
of  the  high  shear  rates  encountered  in  the  pipeline  tests. 
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INTEODUCTION 

The  investigation  into  the  viscosity  characteristics  of  iilberta 
crude  oils  T^-as  initiated  at  the  University  of  Alberta  in  1950.  The 
purpose  of  this  research  program  was  to  obtain  rheological  data 
on  various  Alberta  crude  oils  which  could  be  used  to  predict  the 
power  requirements  of  existing  pipelines  or  used  in  the  design  of 
pipelines  handling  Alberta  ciude  stocks. 

The  research  program  in  chronological  order  has  been: 

(1)  The  laboratory’-  investigations  by  7^.  Sidjak  (17)  on  the 
rheological  properties  of  some  Alberta  crudes  using  the 
MacMichael  viscometer. 

(2)  The  design,  assembly,  and  preliminaiy  testing  of  pilot 
pipeline  using  stove  oil  and  Eedwater  crude  oil  by 

M.  Clmiilesc  (9)  • 

(3)  An  investigation  of  the  viscosity  characteristics  of 
Alberta  crudes  using  the  pilot  pipeline.  The  result 

of  this  investigation  is  presented  in  this  dissertation. 

An  evaluation  of  the  consistency  of  the  liquid  is  necessary 
to  determine  the  pressure  drops  encountered  in  the  pipeline 
flow  of  liquids.  The  consistency  of  a  Nevlonian  liquid  is 
independent  of  the  rate  of  shear  and  is  commonly  known  as  the 
viscosity  of  the  liquid.  For  such  liquids,  pressure  drops  may  be 
calculated  by  standard  foiraulae  from  a  single  evaluation  of 
viscosity  at  the  flo^v  temperature. 


4. 


The  consistency  of  a  non-Ne\7tonian  liquid  is  dependent  on 
the  rate  of  shear,  flo^T  temperature,  and  for  certain  liquids, 
upon  previous  thermal  treatment.  For  these  fluids,  data  relating 
"viscosity”  to  shear  rate  and  tempera. ture  is  needed  for  the 
computation  of  pressure  drops.  Without  such  data,  the  calculation 
of  pressure  drops  becomes  little  better  than  guesswork. 

Most  Alberta  crudes  behave  as  Newtonian  fluids  at  temperatures 
above  the  '’wax  point”.  Previous  investigations  by  Sidjsk  (l7)  and 
Bauer  and  Ruston  (5)  on  some  Alberta  crude  oils  indicate  wax  points 
in  the  region  50  to  60  degrees  F.  ’’Viscosity*'  -  temperature 
relationships  have  be®i  found  to  be  near  linear  above  these 
critical  points.  Below  wa;t  point  temperatures,  however,  some 
crude  oils  tested  have  sho\m  non-Newtonian  characteristics. 

Sidjak  (17)  has  shown  that  these  crudes  at  temperatures  below 
the  '’?7ax  point’’  will  behave  as  pseudoplastic  materials. 
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THEORY  MD  LITERATURE  REVIEW 
Classification  of  Fluids 

Fluids  are  classified  according  to  the  relationship  of 
shear  stress  to  the  rate  of  shear,  i.e.,  the  shape  of  shear 
stress  -  rate  of  shear  curves  as  shovjn  in  Figure  3,  (Alves, 
Boucher  and  Pigford  2) .  A  fluid  subj  ected  to  a  shear  stress 
undergoes  continuous  deformation  and  the  resistance  by  the 
fluid  to  the  applied  stress  is  called  its  consistency. 
Numerically  the  measure  of  consistency  is  the  ra.tio: 

c  gc 

-(dv/dr)  (l) 


where  ^  -  unit  shear  stress, 

-  (dv/dr)  r  rate  of  shear  at  the  point  of  stress. 
Sc  =  conversion  constant 


At  any  given  temperature,  Newtonian  fluids  are  characterized 
by  a  constant  ra.tio  between  the  rate  of  shear  and  the  applied 
shearing  stress  (Figure  3,  curve  A).  For  such  substances,  the 
consistency  is  called  viscosity  and  defined  by  the  equation: 


^  ^  /'civ 


c  - 


(2) 


where  ^  is  the  coefficient  of  viscosity  or  simply, 
viscosity  and  is  indicated  by  the  slope  of  curve  A. 

If  the  consistency  is  dependent  on  the  rate  of  shear  at 
constant  temperature,  the  fluid  is  said  to  be  non-Newtonian. 


I 
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For  these  fluids,  the  property  corresponding  to  the  viscosity 
of  a  Newtonian  has  been  described  in  a  multiplicity  of  terms, 
e.g.  plasticity,  rigidity,  apparent  viscosity  and  plastic 
viscosity.  Non-Newtonian  fluids  have  been  further  classified 
into  distinct  types  such  as  Bingham  plastics,  pseudoplastics, 
thixo tropes,  dilatants,  etc.,  and  present  shear  ra.te  -  shear 
stress  relationships  which  differ  from  each  other  and  from  a 
Newtonian  as  illustrated  in  Figures  3,  4  and  5  (lives  l)  • 

1.  Bingham  Plastic  (Curve  B,  Figure  3) 

The  Bingham  plastic  exhibits  a  yield  stress  T  ,  which  must 
be  exceeded  before  motion  is  imparted  to  the  fluid.  At  values 
greater  than  the  yield  stress,  the  ratio  of  shear  stress  to 
shear  rate  is  constant  and  is  called  the  coefficient  of  rigidity 
or  plastic  viscosity  (Green,  ll) .  The  coefficient  of  rigidity, 

,  is  defined  by  the  eq,UE.tion: 

f  =  ~  %)  Sc  (3) 

-0 

where  -  yield  stress 

Experimental  data  on  Bingham  plastics  give  shear  stress  - 
shear  rate  relationships  similar  to  Figure  4,  snd  at  very  low 
shear  rates,  the  straight  line  relationship  between  shear 
stress  and  shear  rate  degenerates  as  shown. 

Sewage  sludge  and  grain  suspensions  in  water  are  materials 
which  exbibit  Bingham  plastic  beha.viour. 


■*a^. 
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2,  Pseudoplastic  (Curve  C,  Figure  3) 

Pseudoplastic  fluids  flow  under  any  deforming  force,  but 
the  consistency  of  these  materials  decreases  at  increasing- 
rates  of  shear  as  shorn  by  the  decreasing  slope  of  curve  C, 


indicated  by  the  slope  of  the  S'&raight  line,  and  the  slope 
of  the  tangent  at  zero  rate  of  shea.r  indicates  the  viscosity 
at  zero  rate  of  shear,  These  ttro  limits  are  dependent 

on  the  system  and  temperature,  but  are  independent  of  rate 
of  shear  and  the  dimensions  of  the  confining  tube  (Wind±ng, 
Baumann  and  Kranich,  2l) . 

The  consistency  curve  of  a  pseudoplastic  fluid  is  best 
represented  by  a  complex  function: 


(4) 


Jm  approximation  of  the  curve  in  the  region  where  linearity 
exists  between  shear  stress  and  rate  of  shear  is  given  by  the 
equation! 


(5) 


the  effective  coefficient  of  rigidity. 

In  the  evaluation  of  the  consistency  of  pseudoplastic  fluids. 


the  term  "apparent  viscosity"  has  been  used,  and  has  been 
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defined  as  "the  viscosity  of  a  pseudoplastic  if  it  were  a 
Newtonian  material  with  a  consistency  curve  passing  through 
the  point  of  measurement"  (Green,  11).  Bef erring  to  curve 
C  of  Figure  3,  the  apparent  viscosity  of  a  pseudoplastic  is 
indicated  by  the  slope  of  line  0~c  if  c  represents  the  point 
of  measurement. 

Winding,  et  al  (2l)  in  discussing  the  pseudoplastic 
behaviour  of  GRS  latices,  suggest  that  the  suspension  of 
finite  particles  in  a  liquid  may  result  in  decrease  of  apparent 
viscosity  at  high  rates  of  shear.  The  decrease  may  be  attributed 
to  any  one  or  a  combination  of  elonga.tion  of  spherical  particles 
by  disintegration  of  aggregates,  alignment  of  particles  possessing 
rod-like  structures,  or  through  the  breakdowm  of  solvation  layers 
of  fluid  about  the  particle. 

PJhen  temperature  conditions  are  favorable,  precipitation  of 
wax  takes  place  in  many  crude  oils.  These  particles  could 
influence  the  rheological  behavior  of  crude  oils  through  one  or 
more  of  the  mechanisms  previously  suggested.  Sidjak  (17)  found 
that  when  temperature  conditions  resulted  in  wax  precipitation 
the  crude  oil  behaved  in  a  pseudoplastic  manner.  Padgett  (14) 
found  that  petroleum  wax  could  crystallize  out  in  several 
different  foims,  e.g.  as  plate  crystals,  needle  crystals,  or 
in  a  mal-crystalline  form,  end  that  the  tj^pe  of  crystal  formed 
was  dependent  on  the  characteristics  of  the  wrax  present  in  the 
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oil,  the  concentration  of  wax,  the  rate  of  cooling  and  on  the 
influence  of  crystal  forms  upon  one  another.  Investigations  by 
Sidjak  (17)  and  Bauer  and  Raston  (5)  showed  that  apparent 
viscosities  of  certain  Alberta  crades  were  also  influenced  by  the 
theimaJ.  treatment,  to  which  the  oil  has  been  subjected  prior  to 
viscosity  measurement.  The  thermal  treatment  of  the  crude  may  have 
caused  fomation  of  different  crystal  foims. 

3«  Dilatant  Materials  (Curve  D,  Figure  3) 

Dilatant  materials  are  defined  by  the  general  relationship 
between  shear  stress  and  rate  of  shear  as  shown  by  curve  D.  The 
consistency  of  these  fluids  increases  with  increasing  shear  stress. 
A  characteristic  feature  of  the  dilatant  material  is  that  a 
stiffening  of  the  material  takes  place  when  it  is  subjected  to  a 
shearing  stress. 

Dilatancy  is  tj^ical  of  heavy  suspensions  in  which  particles 
have  settled  into  a  position  of  ’'minimum  voids.”  The  disturbance 
of  particles  under  shear  will  expand  or  "dilate”  the  volume  of 
voids  resulting  in  an  insufficiency  of  fluid  vehicle,  consequently, 
a  dry  appearance  and  stiffening  of  the  suspension  takes  place 
(Green,  11) .  A  common  material  which  exhibits  dilatancy  is  wet 
sand  whose  water- to-sand  ratio  is  just  sufficient  to  fill  the  voids 
when  the  voids  are  at  their  minimum  volume. 

Dilatant  behaviour  has  not  been  observed  in  crude  oils. 
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4.  Thixotropic  Materials  (Figure  5) 

Tiiixotropic  materials  possess  a  structure,  the  breakdown  of 
which  is  a  function  of  time  as  well  as  of  the  rate  of  shear.  These 
structures  can  relzuild  themselves,  if  not  prevented  from  doing  so 
by  externally  applied  forces.  As  a  simple  illustration  of 
thixotropy,  the  addition  of  a  suitable  amount  of  electrol^d^e  to 
a  strong  sol  of  ferric  oxide  or  alumina  results  in  a  jelly  if 
allowed  to  stand.  Shaking  the  jelly  will  reform  the  sol  which  will 
set  again  on  standing. 

Data  from  rotational  viscometers  for  thixotropic  materials  give 
shear  stress-rate  of  shear  diagram  similar  to  Figure  5  with  the 
characteristic  hysteresis  loop  which  is  an  indication  of  the 
degree  of  thixotropy  of  the  materiak.  Referring  to  Figure  curve 
abc  is  obtained  wtdle  the  material  is  being  continually  broken 
down,  increased  shear  rate  resulting  in  decreased  consistency. 

If  the  downcurve  ca  is  run  immediately  after  the  completion  of 
upcurve  abc,  no  further  brealcdown  occurs,  h^ce  higher  rates  of 
shear  are  obtained  for  the  corresponding  shear  stresses.  1 
repetition  of  loop  a.bca  is  obtained  if  the  thixotropic  material 
is  retested  after  being  allowed  to  stand  (Alves,  et  al,  2) . 

An  application  of  thixotropic  behaviour  is  foimd  in  drilling 
muds.  Water  alone  is  not  adequate  as  a  drilling  fluid  since  it 
is  not  heavy  enough  and  cuttings  may  easily  settle  out  end 
freeze  the  drill.  Thiro tropic  suspensions  of  bentonite  are  used 
as  boring  fluids  to  maintain  the  required  hydrostatic  head  end 
to  bring  up  the  cuttings.  These  muds  set  to  soft  jels  which 


are  easily  liquidied  (Weiser,19) 
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Other  materials  exhibiting  thixotropic  behaviour  are  paints 
and  printing  inks.  The  thixotropic  behaviour  of  some  cinide  oils 
has  been  reported.  This  is  not  surprising  since  come  crude  oils 
will  gel  if  allowed  to  stand  at  low  temperatures. 


Pipeline  Flow  Equations  for  Newtonian  Fluids 

A.  Eeynolds  Criterion  of  Flow. 

The  flow  of  Newtonian  fluids  may  take  place  in  a  viscous, 
transitional,  or  turbulent  manner.  Viscous  flow  is  characterized 
by  the  movement  of  the  fluid  in  layers  or  laminae,  and  turbulent 
flow  by  the  presence  of  eddying  and  turbulence  in  the  flowing 
stream.  Laminar,  transitional  or  turbulent  flow  is  dependent  on 
the  nmerical  value  of  the  dimensionless  ratio  which  is  called 
the  Reynolds  number  Nj^.  Experimental  results  indicate  that 
laminar  flow  exists  when  values  of  Reynolds  number  are  less  than 
2100  (  Walker,  Lewis,  McAdams  and  Gilliland,  18) .  The  value  of 
2100  is  called  the  lower  critical  value  of  Reynolds  number.  The 
upper  critical  value,  which  indicates  the  lowest  value  of  Reynolds 
number  for  which  complete  turbulence  is  present,  is  reported  by 
Brown  (7)  as  3000  while  Binder  (6)  reports  an  upper  critical  value 
of  4000.  Values  lying  between  these  two  critical  Rejnolds  numbers 
indicate  a  state  of  transition  where  flow  is  neither  laminar  or 
turbulent  but  a  mixture  of  these  txTo  states. 
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B.  The  Flow  Equation 

For  the  solution  of  problems  involving  steady  state  flow  of 
essentially  non compressible  fluids.  Walker  et  al  (l8)  present  the 
mechanical  energy  equation  in  the  foim: 

-  /f®  =  ‘^2  -  k  t  Vj-Vf  +  P  (6) 

2ge 

where  =  external  work  done  on  the  flowing  fluid 

0  =  specific  volume  of  flowing  fluid 

X  =  height  above  datum. 

P  =  pressure 

V  =  velocity 

g  -  acceleration  due  to  gravity 
S’c  ~  conversion  constant 
F  =  total  friction  due  to  flow. 

In  the  case  of  fluid  flow  in  a  horizontal  pipe  of  constant  cross- 
section  where  change  in  specific  volume  of  the  fluid  is  very  small, 
and  in  the  absence  of  external  work,  equation  (6)  reduces  to: 

-AZ  -  F  (7) 

Pav 

where  -  average  density  of  fluid. 

The  friction  tem  in  equation  (7)  may  be  evaluated  by  the 
Fanning  equation  which  expresses  pressure  drop  in  terms  of  fluid 
density,  flow  rate,  pipeline  geometry  and.  a  friction  factor; 
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e  2gcD 

where  D  =  internal  diameter  of  pipe 
L  =  length  of  pipe  section 
f  -  friction  factor 

In  seme  engineering  texts  equation  (8)  is  presented  as 
^  P  -  2f“  v^L 

p  - — 

gcD 

j  where  f’  =  ^ 

1  ^ 

I 

I  In  this  dissertation  the  values  of  friction  factor  will  be  defined 

I 

I  by  equation  (8). 

C.  Friction  Factor  Equations 

Experimental  data  ha.ve  shown  that  the  friction  factor  is 
dependent  on  two  dimensionless  ratios,  the  RejTiolds  number^ 

;  (Wp)  ^  and  the  relative  roughness ,  (  e  ) «  i.e,,  the  ratio  of  the 

absolute  roughness,  (e),  to  the  diameter  of  the  pipe,  (D). 

In  laminar  flow,  relative  roughness  of  the  pipe  wall  has  no 
;  effect  upon  the  friction  factor  and.  the  experimentally  determined 

'  relationship  between  Reynolds  number  and  friction  factor  is: 

^  ^  (9) 

“e 

where  -  DVp«  Reynolds  niunber 

;  Poiseuille's  equation  states  that  the  pressure  drop,  encountered 

I  in  the  laminar  flow  of  a  fluid  through  a  tube,  is  directly 

proportional,  to  the  product  of  viscosity,  average  velocity  and 
length  and  inversely  proportional  to  the  square  of  the  diameter. 


15 


that  is: 


This  equation,  while  derived  theoretically,  is  equivalent  to 
the  combination  of  equations  (S)  and  (9). 

Empirical  foimulae  have  been  developed  to  show  the  relationship 
between  friction  factor,  Reynolds  number  end  relative  roughness  in 
turbulent  flow.  For  smooth  tubes  the  equations  are  given  in  terms 
of  Reynolds  number  and  friction  factor.  Equations  commonly  used 
for  smooth  tubes  are: 

Blasius  equation  (4) 


f  =  0,316 


(10) 


Koo  equation  (15>  18) 

f  =  ,0056  0,50 


(11) 


Kaiman,  Nikuradse  and  Prandtl  equation  (6,l6) 


/f 


3^  =  2  log^Q  \ 


1 


(12) 


Stanton  and  Pannell  equation  (3) 

f  =  0,3,'"05 

(NgV-252 


(13) 


Heltzel’s  equations  (3,12) 


f  =  0,364 


(14) 
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^  =  0*157 

N  0.188 


(15) 


R 


The  Blasius  equation  is  especially  suited  for  Reynolds  numbers 
ranging  from  10,000  to  100,000.  The  Koo  equation  is  based  on  a 
Reynolds  number  range  from  3000  to  3,000,000.  Anderson  (3)  showed 
that  the  Stanton  and  Pannell  relationship  could  be  used  to  compute 
friction  factors  in  the  Reynolds  niamber  range  of  6000  to  130,000 
for  the  transmission  of  some  crude  oils  through  large  diameter 
lines  (l6,20  and  30  inch).  Equations  (14)  and  (l5)  were  the  results 
of  Heltzel's  work  on  lines  6  to  12  inches  in  diameter.  He  found 
that  friction  factors,  calculated  in  the  flow  of  crude  oil  through 
these  lines,  could  be  closely  represented  by  equation  (14)  for 
values  of  Reynolds  numbers  in  the  range  2,500  to  57,600  and  by 
equation  (15)  for  values  of  Reynolds  numbers  greater  than  57,600. 

In  rough  pipes,  friction  factors  encountered  in  turbulent  flow 
are  influenced  by  the  value  of  relative  roughness  as  well  as  the 
Reynolds  number.  Furthermore,  the  friction  factor  approaches  a 
limit  dependent  on  the  relative  roughness.  Karman's  equation  for 
the  value  of  this  limit  in  fully  developed  turbulent  flow  is 
given  by  Rouse  (l6)  as: 


(16) 


In  general,  pressure  drops  in  fluid  flow  are  computed  by  using 


the  Fanning  friction  factor  equation  and  by  evaluating  friction 
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factors  from  the  standard  friction  factor  -  Reynolds  number  plots 
given  in  fundamental  engineering  texts.  Such  plots  present  in 
graphical  form,  equation  (9)  for  laminar  flow  and  the  equivalent 
of  (12)  for  turbulent  flow  in  smooth  pipes.  Parameters  indicating 
lines  of  constant  relative  roughness  show  the  effect  of  the  wall 
roughness  on  the  friction  factor. 

D.  Shear  Stress  and  Shear  Rate 

r 

>  a 


Wall 


FI6.  6  5Kear  Stress  Distribution  in  Circular  Conduit 

In  the  isothermal  steady  state  flow  of  a  cylindrical  mass  of 
fluid,  consider  the  section  of  radius  r  and  pressure  drop, 

=iJP  (Figure  6).  A  force  balance  may  be  written: 

=  (27,  rL)t'  (17) 
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'|r  ^P)  represents  the  difference  in  force  between 

the  two  faces,  and  2  ii  rLC represents  the  opposi-ng  force  acting  on  the 
cylindrical  surface  of  radius  r.  The  shear  stress  is  shown  to  be: 

-  rAP  (18) 

2L 

Similarly  the  shear  stress  at  the  wfill  may  be  shown  to  be: 

tw  =  D— (^-  (19) 

4L 

Since  no  limitations  were  made  as  to  the  rheology  of  the  fluid 
or  the  type  of  flow  (laminar,  transitional  or  turbulent),  equation 
(l8^  shows  that  the  steady  state  flow  of  all  "fluids  will  be  in  such 
a  manner,  that  a  linear  shear  stress  distribution  will  be  obtained 
(Figure  6). 

For  Newtonian  fluids  in  laminar  motion^  the  shearing  force  at  the 
wall  of  the  pipe  is: 


Substituting  this  in  equation  (19),  the  rate  of  shear  at  the  iwll 
for  a  Newtonian  fluid  in  laminar  motion  is: 

(20) 

Fluid  flow  through  tubes  and  pipes  is  characterized  by  a 
velocity  profile  with  maximum  velocity  at  the  axis  and  zero 
velocity  at  the  wall.  A  boundary  layer  exists  adjacent  to  the 
Awall  in  which  fluid  is  flowing  in  a  laminar  fashion  while  the  main 
core  of  the  stream  may  be  turbulent,  (Coulson,  10).  The  shear 


*#1 
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A-  U 
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rate  at  the  wall  is  again  given  by  equation  (20)  because  of  the 
laminar  boundary  layer. 


Pipeline  Flow  Equations  for  Non-Newtonian  Fluids 

The  mechanical  energy  equation  (6)  and  the  Fanning  equation  (8) 
may  be  used  to  estimate  the  pressure  drops  encountered  in  the 
pipeline  flow  of  non-Newtonian  fluids  as  well  as  Newtonian  fluids. 

However,  for  non-Newtonians  the  equations  relating  the  friction  factor 
to  the  Reynolds  number  and  the  relative  roughness  for  Ne?rbonian  fluids 
are  no  longer  valid.  If  an  apparent  Reynolds  number  is  defined  for 
the  flow  of  non-Newtonians,  using  some  other  consistency  term  for 
viscosity; 

apparent  N^^  =  YNp  (21) 

(consistency  of  non-Nevrbonian) 

it  is  found  that  the  apparent  Reynolds  number  will  vary  throughout 
the  cross  section  of  the  pipe  since  the  consistency  of  non-Newtonians 
are  dependent  on  the  rate  of  shear  which  varies  from  maximum  at  the 
pipe  wall  to  zero  at  the  axis  of  the  pipe.  Establishment  of  the 
character  of  flow,  i.e.  laminar,  transition  or  turb\ilent,  has  been 
done  by  most  investigators  by  substituting  coefficient  of  rigidity, 
limiting  viscosity  or  apparent  viscosity  for  the  consistency  term  in 
the  apparent  Reynolds  number. 

Buckingham  (ll)  in  1921  derived  a  theoretical  expression  for  the 

A 


capillary  flow  of  a  plastic  based  on  the  definition  of  a  plastic. 
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equation  (3) •  He  pictured  "plastic  flow"  to  be  made  up  of  three 

parts: 

(1)  the  flow  of  a  central  plug  in  which  no  shear 
occurred. 

(2)  a  laminar  layer  enveloping  the  plug  and 

(3)  a  lubricat-i-ng  layer  adjacent  to  the  wall  with  a 
finite  slippage  velocity  v^^  and  thickness,  £-  . 

Buckingham  assumed  that  the  thin  lubricating  layer  e:>iiibited  a 
Newtonian  viscosity.  The  final  form  of  his  equation  is: 

Q  =  (ap  -  4ap  i  (22) 

8/jL  3  3^P^  2f^L 

volume  rate  of  flow 
radius  of  capillary 
coefficient  of  rigidity 

length  of  capillary 
total  pressure  drop 
pressure  required  to  initiate  motion 
viscosity  of  a  Nevrtonian  fluid 

thickness  of  lubricating  layer 
Equation  (22)  shows  that  irhen  4P  =  isp,  the  first  term  becomes 
zero,  and  the  flow  becomes  entirely  due  to  slippage.  Whennp  is 
less  than  n^p , there  is  no  flow  since  the  total  applied  stress  is 
less  than  the  yield  stress.  When  there  is  no  slippage  the  last 
tem  vanishes  and  plug  flow  is  not  eliminated  completely  until 
the  volume  rate  of  flow  becomes  infinite.  Although  Buckingham’s 


where  Q  - 
R  = 

^  = 
L  = 

AP  = 

AP  = 

t  = 
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theoretical  equation  is  important  in  discussing  the  plug  and 
laminar  flow  of  plastics,  the  equation  is  not  satisfactorily 
solved  due  to  the  numbrous  unknown  quantities  inherent  in  it.  A 
simpler  flow  equation  by  Bingham,  again  based  on  the  definition  of  a 
plastic  and  cited  by  Babbitt  and  Ga-ldxvell,  (8)  is; 

Q  =  ITR^go  (^'P  -  Aip  +  1  V-,)  (23) 

■^L  3  3 

Bingham  did  not  assume  slippage  at  the  wall  but  only  assumed  a 
central  plug  moving  inside  a  laminar  layer.  In  terms  of  velocity, 
the  equation  becomes: 


V  = 


°c 

32y^L 


(^P  -  A-p  X  1 

3  3 


(24) 


Substituting  equation  (24)  in  the  Fanning  equation  results  in  the 
friction  factor  expression: 

f  (2gcD)f2M)2(aP  -  I  1  (25) 

(°  L  Wga’  3  3 


A  dimensional  analysis  by  Winning  (22)  on  the  properties 
assumed  to  influence  the  pipeline  flow  of  plastics  resulted  in  the 
expression: 


'  f  D 


(26) 


From  Winning'  s  experimental  data,  a  plot  of  INp  vs.  friction  factor 
showed  the  effect  of  the  "yield  ratio",  D  ^  g  ,  in  the  laminar 
flow  region.  For  values  of  JNp  to  2000,  lines  of  constant  yield 

/f 

ratio  were  obtained  parallel  to  the  equation  f  =  64  (z^) . 

iNp 


In  this 
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region,  higher  va].ves  of  yield  ratio  corresponded  to  higher  values 
of  friction  factor  for  any  given  value  of  W/f  and  the  line 


f  =  64  )  represented  a  yield  ratio  value  of  zero. 

DV/? 

Another  method  of  evaluating  friction  factors  in  the  laminar  flow 
of  non-Ne?rbonian  fluids  has  been  through  the  use  of  ’’apparent 
viscosities”.  A  modified  Poiseuille  equation  in  which  apparent 
viscosity  has  been  substituted  for  viscosity  is: 


aP  = 


LV 


(27) 


where^^  is  the  apparent  viscosity  and  depicts  the  net 
result  of  the  consistency  properties  of  the  non-Newtonian  in  laminar 
flow  through  a  tube.  Substitution  of  equation  (27)  in  the  Fanning 
equation  (8)  gives; 

f  =  64  (^SL  )  (28) 

cy 

or  f  =  64 

(«E), 

where  (N^^  is  apparent  Reynolds  number. 

A 

In  fully  developed  turbulent  flow,  theoretical  equations 
describing  the  flow  of  non-Newtonian  fluids  ar*e  not  Imown. 
Experimental  results  of  Winning,  (22),  Alves,  Boucher  and  Pigford  (2) 
Winding,  Baumann  and  Kranich  (21)  show  that  non-NeTrbonians  beha.ve 
in  a  manner  similar  to  Newtonians  in  the  turbulent  region,  that  is. 


the  effective  consistency  is  not  influenced  by  the  shear  rates  of 
the  non-Nevrbonian  fluids.  The  work  of  Alves  et  al  (2)  on  the  turbulent 
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flow  of  a  plastic,  an  eighteen  percent  titanium  dioxide  slurry, 
over  a  range  of  apparent  Reynolds  numbers  of  4,000  to  100,000 
gave  values  of  apparent  viscosities  of  1.5  to  1.6  centipoise. 

With  these  results  they  showed  that  the  conventional  friction 
factor  -  Reynolds  number  charts  could  be  used  to  obtain  values 
of  friction  factor  using  apparent  Reynolds  numbers.  Winding, 
Baumann  and  Kranich  (21)  suggest  tha,t  limiting  viscosity  at 
infinite  shear,  in  the  case  of  pseudoplastics  or  the  coefficient 
of  rigidity,  in  the  case  of  plastics,  can  be  used  to  compute 
apparent  Reynolds  numbers  in  turbulent  flow  and  that  these 
apparent  Reynolds  numbers  can  be  used  to  determine  friction 
factors  from  the  standard  friction  factor  -  Reynolds  number 
charts. 

The  literature  review  indicates  tha-|i  in  fully  developed 
turbulent  flow  of  non-Nevrtonian  fluids,  friction  factors  may 
be  evaluated  by  the  method  of  using  apparent  Reynolds  numbers. 

A  suitable  friction  factor  equation  such  as  the  Karman  Prandtl 
equation  (l2)  or  the  Heltzel  equation  (l4)  may  be  used  in  relating 
the  apparent  Reynolds  number  to  the  friction  factor.  For  example, 
modification  of  the  Heltzel  equation  becomes: 


(29) 
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As  in  the  modified  Poiseuille  equation,  the  apparent  viscosity 
terra  in  equation  (29)  depicts  tlie  net  effect  of  the  consistency 
property  of  the  non-Ne\7tonian  fluid  in  fully  developed  turbulent 
flow. 

A  few  friction  factor  equations  have  been  presented  for 
non-Newtonian  fluids.  Since  Alberta  crude  oils  have  been  found 
to  behave  in  a  pseudoplastic  manner  at  low  temperatures  (Sidjak,  17)^ 
the  Buckingham,  Bingham  and  Winning  equations  for  plastic  flow 
cannot  be  used.  The  method  of  apparent  viscosities  seems  the 
effective  method  for  evaluating  friction  factors,  especially 
for  the  case  of  the  pilot  pipeline. 

N  omen  clature 

Definition  of  Symbols  Used  in  Theoiy  and  Literature  Review. 

-  ^  -  velocity  gradient  or  rate  of  shear,  sec 

dr 

-  rate  of  shear  at  the  pipe  wall,  sec 

^dr^w 

D  -  internal  diaiiieter  of  pipe,  ft. 
e  -  absolute  roughness  of  pipe  wall,  ft. 
f  =  fanning  friction  factor,  dimensionless. 

F  -  friction  loss,  3.b^  -  ft. 
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g  z  acceleration  due  to  gravity,  ft/ sec^ 

g^,  z  conversion  constant  =:  32.2  -  ft. 

Ib^  -  sec^ 

L  =  length  of  pipe,  ft. 

z  Reynolds  number,  dimensionless 


=  apparent  Reynolds  number,  dimensionless 
a 

_  2 
P  z  pressure,  lb/  ft. 

aP  =  pressure  drop  required  to  initiate  flow, for  plastic^  Ib/ft, 
z  total  pressure  drop,  Ib/ft. 

Q  r  volume  rate  of  flow,  ft^/ sec, 
r  =  radial  distance  from  pipe  axis,  ft. 

R  -  internal  radius  of  pipe,  ft, 

V  -  velocity  of  fluid  ft/sec. 

V  -  point  velocity  of  fluid,  ft/sec, 

Vj^  z  slippage  velocity  of  plastic  fluid,  ft/sec, 
w  =  external  work  done  on  floirilng  fluid,  lb£>  -  ft. 

X  z  height  above  datum,  ft. 


lb 


m 


thickness  of  lubricating  layer  in  the  Buckingham  equation, 
function 

viscosity,  Ib/ft-sec. 
apparent  viscosity,  Ib/ft-sec. 


ft. 
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-  limiting  viscosity  of  pseudoplastic  at  zero  shear,  Ib/ft-sec. 
z  limiting  viscosity  of  pseudoplastic  at  infinite  shear,  Ib/ft-sec. 
=  coefficient  of  rigidity,  Ib/ft-sec. 
z  effective  coefficient  of  rigidity,  Ib/ft-sec* 

-  fluid  density,  Ib/ft.^ 

3  shear  stress,  Ib/ft.^ 

shear  stress  at  the  wall,  Ib/ft.^ 

=  yield  stress  for  plastic,  Ib/ft. 

=  effective  yield  stress,  Ib/ft.^ 

l)  -  specific  volume  of  fluid,  ft^/lb. 
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DESCRIPTION  OF  EXPERIMENTAL  EQUIPMENT 

The  essential  features  of  the  pilot  pipeline  (Figures  1 
■aid  2)  are  the  refrigerated  shelly  ref rigeration  unit, test  lines, 
precooling  lines,  crude  oil  storage  and  measurement  tanks, 

Roper  gear  pump,  and  the  control  instruments. 

T\to  lengths  of  eight  inch  pipe,  fitted  mth  dead  end 
flanges  (Figure  7),  maie  up  the  outer  casing  of  the  refrigerated 
shell.  An  eight  foot  length  encloses  the  calming  section  and 
an  eighteen  and  a  ha.lf  foot  length  encloses  the  test  section. 

The  casings  a.re  connected  by  vapor  and  liquid  eq-uilibrium  lines 
(Figure  S) .  The  entire  casing  is  insulated  Tsdth  glass  -cvool  and 
covered  r/ith  a  vapor  barrier  of  corrugated  aluminum  sheet  to 
prevent  heat  transfer  from  the  surroundings  into  the  shell. 
Arrangement  of  test,  refrigeration  and  precooling  lines  Tlthin 
the  shell  is  shofmi  schematically  in  Figure  2.  Test  and  precooling 
lines  are  immersed  in  a  bath  of  liquid  Freon-12.  A  pressuretrol, 
■which  is  actuated  by  the  vapor  pressure  of  Freon-12  -within  the 
shell,  controls  the  bath  temperature.  Settings  on  this  controller 
turn  on  strip  heaters  if  the  pressure  decreases  and  shut  off  the 
heaters  if  the  pressure  increases.  By  this  means,  a  relatively 
constant  pressure,  and  hence  constant  temperature  may  be  maintained 
within  the  shell.  Other  accessories  to  the  shell  include  the 
liquid  level  gauge,  pressure  gauge,  therraometer  and  thermometer 


well 


I 


28 


Three  stainless  steel  test  lines  are  mounted  on  the  unit 
as  sho\m  in  Figure  7.  The  internal  diameters  of  these  lines  are 
0.249  inches,  0.311  inches  and  0.386  inches.  The  manifold 
arrangement  for  selection  of  the  test  line  is  shorai  in  Figures 
10  and  11.  The  test  section  beteeen  pressure  taps  is  19.5  feet 
and  is  preceded  by  a  calming  section  of  8.5  feet,  Weston  dial 
thermometers,  of  range-40  to  140 °F.,  are  inserted  in  each  of  the 
distributing  manifolds  to  measure  the  flow  temp  era.  ture  of  the 
crude  oil.  Figure  13  shows  a.  tj^ical  pressure  tap  assembly. 

Care  was  taken  in  fabrication  to  remove  any  burrs  which  could 
induce  turbulence  at  the  pressure  taps. 

Three  pressure  gauges  are  used  to  deteimiine  the  pressure 
dxop  over  the  test  section j  an  upstream  Heise  gauge  of  range 
0  to  1000  psi.  in  2  psi.  intervals,  a  do^mstream  Acragauge  of 
range  0  to  100  psi.  in  one-half  psi.  intervaJ-s,  and  a  mercusry- 
under-oil  differential  manometer  for  lov;  pressure  drops  (O  to  30 
psi. ). 

Two  60  gal.  capacity  tanks  (Figure  12),  mth  sight  glasses 
provide  for  the  volumetric  measurement  of  the  fluid  used  in  the 
tests.  The  valves  and  pipes  are  arranged  so  that  either  tank 
may  be  used  as  the  reservoir  or  the  measuring  tank.  Drain.-,  ge 
pipes  are  installed  for  manual  draining  of  tanks  of  fluid  which 
cannot  be  pumped  out.  The  tanks  are  completely  enclosed  to  prevent 
weathering  of  the  crude,  however  a  safety  relief  valve  is  used 
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to  prevent  build-up  of  excessive  tank  pressures, 

A  Roper  gear  pump,  capable  of  delivering  pressures  up  to 
1000  psi.  is  used  to  circulate  the  crude  oil  through  the  system. 
Piping  and  valves  are  arranged  to  permit  part  or  all  of  the  crude 
oil  to  be  returned  to  the  storage  tank  before  or  after  precooling 
of  the  oil  (Figure  2)  •  l?hen  part  of  the  oil  stream  is  by-passed 
before  cooling,  the  temperature  of  the  storage  tank  oil  remains 
essentially  at  the  higher  temperatures  (50  to  70  F)  since  any 
cooling  of  the  storage  oil  is  due  to  the  return  of  cold  test  oil 
to  storage.  TiJhen  pai*t  of  the  oil  stream  is  by-passed  after 
precooling,  greater  cooling  of  the  storage  tank  temperature 

Q 

is  effected  and  tank  temperatures  may  reach  as  low  as  35  F, 

A  three  horsepower  Frecn-12  vapor  compression  unit  is  used 
to  refrigerate  the  shell  (Figure  9) .  Freon-12  vapor  taken  on  the 
suction  side  is  compressed,  cooled,  and  liquified,  and  expanded  from 
the  half  inch  liquid  line  into  a  l-g-  inch  copper  refrigeration 
line  which  passes  through  the  length  of  the  shell.  Cooling  of 
the  bath  is  accomplished  by  the  vapors  condensing  on  the  refrigeration 
lines.  This  results  in  the  lowering  of  vapor  pressure  and  temperature 
of  the  liquid  Freon  in  the  bath.  The  crude  oil  is  cooled  by  the 
transfer  of  heat  through  the  walls  of  the  precool  and  test  lines 
into  the  liquid  Freon  bath. 

Pressure  gauges,  magnetic  controls  c-ond  manual  switches  for 
the  heaters,  pump  and  refrigerating  unit  are  mounted  on  an  instrument 


panel  board  as  shown  in  Figure  14. 
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Figure  8.  Equilibrium  Lines  and  Upstream  Pressure  Tap 


Figure  9.  View  of  Upstream  Portion  of  Pilot  Pipeline 

Showing  Refrigeration  Unit  and  Calming  Section 
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Figure  10.  End  View  of  Upstream  Manifold 
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Figure  11.  End  View  of  Downstream  Manifold  Figure  12.  View  of  Crude  Oil 

and  Pressure  Taps  Measuring  Tanks,  Freon 

Storage  Tanks  and 
Instrument  Panel 
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Figure  14.  Instrument  Panel  for  pilot  Pipeline 
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OPERilTION  OF  TI-IE  PILOT  PIPELINE 

The  pilot  pipeline  was  drained  of  all  oil  used  in  previous  tests 
;  before  it  charged  with  new  crude  oil.  Draining  was  followed  by 

two  successive  flushes.  In  the  first  flush,  15  gal.  of  clean  stove 

1 

oil  was  circulated  through  the  system  at  high  velocities  in  order  to 
loosen  and  remove  any  v/ax  which  may  have  been  deposited  on  the  walls 
of  the  test  lines.  This  was  followed  by  a  flushing  with  15  gal.  of 
the  new  crude  oil  to  be  tested. 

In  preparation  for  testing,  approximately  £0  gals,  of  the  crude 
to  be  tested  were  charged  in  the  storage  tanks.  The  crude  oil  was 
then  circulated  tlirough  the  entire  unit  for  at  least  one-half  hour 
to  mix  any  residual  oil  from  the  flushing  ivith  the  charged  crude. 

Since  crude  oil  density  data  were  necessary  for  computations,  a. 
sample  of  the  crude  in  the  unit  was  taken  and  densities  in  the 
temperature  range  of  20  to  75‘^F.  were  determined  by  a  chain  gravitometer. 

The  pilot  pipeline,  when  not  in  use,  wa.s  kept  refrige3?ated  in  order 
to  minimize  loss  of  Freon-12  in  the  bath  th-rcugh  small  leaks  in  the 
I  shell.  For  the  preliminary  operations  of  flushing,  charging  end 

circulating,  the  refrigeration  unit  was  cut  out  to  enable  the  system  to 
;  attain  equilibrium  at  approximately  70°f, 

For  the  first  test  on  the  new  crude,  the  storage  tank  (Figure  2) 
was  charged  vrith  the  bulk  of  the  crud.e  oil  and  the  by-pass  valves 
j  adjusted  to  presj^ure  drops  which  gave  the  desired  rates  of  flow. 

i'  Circulation  of  oil  was  maintained  until  an  oil  temperature  of  75  F. 

; 

_ 


was  indicated.  At  this  temperature,  the  oil  circulation  was  stopped 
momentarily  and  the  tank  valves  adjusted  for  the  measurement  of  test 
oil  in  the  receiver  tank.  The  time  required  to  collect  a  pre¬ 
determined  volume  of  crude  oil  was  recorded.  While  collecting  this 
volume  of  oil,  upstream  and  downstream  pressures,  upstream  and  down¬ 
stream  temperatures,  temperatures  of  storage  oil,  oil  at  the  pump 
discharge  and  the  final  temperature  of  the  crude  in  the  receiver 
were  recorded.  About  three  readings  of  each  of  the  aforementioned 
variables  were  taJcen  during  any  s:ingle  test  to  obtain  average 
temperatures  and  pressures.  After  the  test,  the  measured  crude  was 
transferred  back  into  the  storage  tank. 

For  the  second  test,  the  refrigeration  unit  was  put  into  operation 
and  the  oil  circulated  in  the  unit  until  the  oil  was  cooled  to  70°F, 
A.t  70°F,,  the  procedure  of  the  first  test  wus  repeated.  Successive 
tests  were  made  at  each  5^F,  drop  in  temperature  to  a  low  of  about 
20®F. 

For  test  temperatures  down  to  35  or  40°F,,  the  oil  was  by-passed 
after  precooling,  this  resulted  in  cooling  of  the  storage  oil  to 
about  40*^F.  It  was  found  that  35°F.  mus  near  the  minimum  attainable 
test  temperature  with  this  arrangement,  even  at  very  low  rates  of 
flow. 

In  acquiring  data  at  lower  temperatures,  the  oil  was  by-passed 
before  precooling.  The  resulting  Ioiy  rates  of  flow  in  the  precool 
lines  made  possible  minimum  temperatures  that  ivere  from  15  to  20°f, 
lower  than  by  the  other  arrangement. 
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RESULTS 

Pilot  pipeline  tests  were  run  on  the  following  samDles; 

1.  Redwater  crude,  sample  no.  1. 

2.  Redwater  crude,  sample  no.  2. 

3.  Stettler  D2-D3  blend. 

4.  Duhamel  -  Malmo  -  New  Norway  blend. 

5.  Leduc  Woodbend  blend. 

6.  Nisku  blend. 

Information  on  the  source  of  samples,  date  received,  A.P.I. 
gravity,  and  experimental  tests  on  the  crude  oils  is  listed  in 
Table  IIB  of  Appendix  B.  In  Table  IIIB,  the  experimental  density 
data  for  all  crude  oil  samples  in  the  temperature  range  20  to 
75°F.  is  recorded.  This  data  is  pressnted  in  Figure  15  as  a  plot 
of  density  vs.  temperature  of  the  crude  oil. 

The  basic  experimentaJL  data  in  Appendix  B  were  obtained  from 
pilot  pipeline  tests.  The  variables  measured  were: 

1.  temperature  of  storage  tank  oil 

2.  temperature  of  oil  at  pump  outlet 

3.  temperature  of  receiver  tank  oil 

4.  temperature  of  Freon-12  bath 

5.  temperature  of  oil  at  test  inlet 

6.  temperature  of  oil  at  test  outlet 

7.  upstream  pressure 


density  of  stock  in  lbs  /ft. 
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8.  do^mstream  pressure 

9.  volume  of  oil  collected  in  test 

10.  time  for  collection  of  measured  volume  of  oil. 

Tables  IV-B  to  IX-B  show  the  tabulation  of  the  pilot  pipeline 


data  with  the  following  modifications; 

(a)  Temperature  of  oil  at  test  inlet  and  temperature  of 
oil  at  test  outlet  was  recorded  as  average  test 
temperature.  The  average  temperature  was  computed  by 
using  equation  (6-a)  of  Appendix  A. 


i.e.  T 


av 


Tout  -  £££ 
59.2 


^'^out 


(b)  The  pressure  difference  between  upstream  pressure  and 
downstream  pressure  rather  than  the  two  pressure 
readings  was  recorded. 

(c)  Instead  of  volume  collected  and  time  of  collection, 
the  flow  rate  in  cubic  feet  per  minute  was  reported. 
These  flow  rates  were  corrected  to  flowing  tempera.ture 
as  shown  by  sample  calculations  in  Section  1  of 
Appendix  A. 

Appendix  C  presents  the  detailed  results  of  computations  for  the 
experimental  data.  Appendix  B.  Table  I-C  gives  a  concise  summary 
of  formulae  used  in  computations  of  pipeline  data.  These  equations 


are  derived  in  Section  2  of  Appendix  A  and  their  use  is  shown  in  the 
sample  computations  of  Section  3>  Appendix  A.  Tables  II-C  to  VII-C 
present  the  detailed  results  of  computations  for  the  crude  oils  tested 


I 


I 


■  ) 


e 


‘I 

i; 

i 

i 


ii 


* 


ii 


42 


by  pilot  pipeline  and  the  quantities  listed  are: 

1.  average  temperature  of  crude  oil  in  test 

2.  fluid  density 

3.  fluid  velocity 

4.  friction  factor 
5-  Reynolds  number 

6.  apparent  viscosity 

7.  shear  rate  at  Trail 

8.  remarks 

The  fluid  density  was  obtained  from  Figure  15  for  the  particular 
crude  oil  at  the  average  temperature  of  test. 

Fluid  velocity  is  given  by  equation  (8-a)  of  Appendix  A  as: 


X  10^  ft/sec,  for  the  0.249  inch  line 

2.03 


Velocities  for  the  0.311  and  0.386  inch  lines  were  computed  by: 


^  X  10^  ft/ 

3.17  ' 


sec.  for  the  0.311  inch  line 


“  0 _  X  10^  ft/sec.  for  the  0.366  inch  line 


4.87 

The  friction  factors  were  computed  for  fluid  flow  through  the 
0.249  inch  pipe  by  equation  (9-a)  of  Appendix  A: 
f  =  988 (-ap)  ’ 


Equivalent  friction  factor  equations  for  the  other  lines  are: 


f  =  12.35(^P)» 


(0,311  inch  line) 
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f  =  15»32  (^P)« 


(0.386  inch  line) 
where  aP*  is  the  pressure  drop  in  Ibs/in^ 


Apparent  Reynolds  numbers  were  computed  as  shown  in  Section  3  of 

Appendix  A  from  the  equations: 

(H  ^  =  &  for  laminar  flow 

^  a  f 


and  ( H^)  ( 0,. 3 64) ^ turbulent  flow. 


The  following  equations  were  used  to  calcula.te  the  apparent  shear 
rate  at  the  wall,  fdvl ; 

lamina.r  flow 


9.74  Q*  X  10-^  gee 

5.10  Q»  X  10^  sec 


r 

turbulent  flow 

1.231(^')  sec'll 

(0.249  inch  line) 

-1 

1.541(^  sec"^ 

(0.311  inch  line) 

1.912('^')  sec"^ 

(0.386  inch  line) 

■o 

where  Q'  is  in  ftVmin 
and  aP'  is  in  Ib/in.^ 


Apparent  viscosities, calculated  by  the  use  of  the 

follo5¥ing  equations; 

laminar  flow  turbulent  flow 

-3  ..  n  non®  lb  (0.249  inch  line) 

(%l  ft. -sec. 

(0.311  inch  line) 


3.21  X  10“-^  (Ap«)  lb  0.0208 

V  ft.  sec. 


4.99  X  10-3  lb  0.0259  (V)-  lb 

V  ft.  sec.  f't^sec. 

7.71  X  10-3  (^P')  ib  0.0322  lb 

V  ft.  sec.  ft  .sec. 


(0.386  inch  line) 


i 
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The  turbulent  or  laminar  flow  of  fluid  was  noted  under  the 
column  ’’Remarks”.  Wlaere  uncertainty  of  flow  was  present,  computations 
were  made  for  both  laminar  and  turbulent  flow  and  entered  in  the 
tables. 

Table  VIII-C  shows  a  comparison  of  the  apparent  viscosities  obtained 
for  Redwater  samples  no.  1  and  no.  2  from  the  smoothed  plots  of 
Figures  17,  IS,  and  19.  For  sample  no.  1,  apparent  viscosities  in 
the  range  20  to  S0°p^  shown  for  the  0,249  inch  line,  the  0.3S6 
inch  line  and  the  average  of  the  two  lines.  The  deviation  of  each 
from  the  avera.ge  curve  is  given  as  a  percentage.  The  apparent 
viscosity  of  Redxmter  sample  no,  2  is  also  compared  against  the 
average  apparent  viscosities  of  sample  no,  1,  and  the  percent 
deviation  listed  in  the  last  column. 

Table  IX-C  presents  smoothed  apparent  viscosity  temperature 
correlation  for  all  the  oils  tested.  These  values  were  obtained 
from  Figures  17,  18,  19,  20,  21,  22  and  23* 

Table  X-C  shows  the  conversion  of  apparent  viscosity  to 
kinematic  and  Saybolt  Universal  viscosities.  Sample  computations 
used  in  these  conversions  are  given  in  Section  4  of  Appendix  A. 
Equations  for  the  conversion  of  apparent  viscosities  are: 

Equation  (iS-a)  of  Appendix  A  :  |^- 

Equation  (21-a)  of  Appendix  A  :  r  6,72^  (0.0022©)  -  1.8  ) 

^  9 


where  )<  =  kinematic  viscosity  in  ft  /sec. 
0  Saybolt  universal  seconds. 
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Appendix  D  shows  the  refrigeration  studies  conducted  on 
the  pilot  pipeline  and  the  context  of  the  Appendix,  is; 

1.  a  brief  introduction  of  the  purpose  of  this  investigation 

2.  experimental  data  for  refrigeration  studies  (Table  I-D) 

3.  detailed  outline  of  computations 

4.  presentation  of  results  (Table  II-D) 

5.  discussion  of  results. 
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DISCUSSION  OF  RESULTS 

A.  Thermal  Behaviour  of  Crude  Oil  in  Test 

Tests  conducted  on  the  pilot  pipeline  show  that  crude  oil 
is  subject  to  distinct  temperature  variations,  depending  on 
the  particular  flow  circulation  selected  for  testing.  Test 
no,  l6  on  the  Duhamel-Maljip.o-New  Norway  blend  was  selected  to 
illustrate  the  temperature  variations  encountered  during  a 
series  of  tests,  since  a  ^/ide  variation  of  test  temperatures 
was  obtained.  The  thermal  behaviour  of  Duhamel-MaLmo-New 
Norv^ay  blend  in  test  no.  16  is  shown  in  Figure  16. 

For  case  I,  where  all  the  oil  pumped,  was  sent  through 
precool  and  then  by-passed  before  test,  the  storage  temperature 
decreased  in  a  regular  fashion  for  each  successive  run.  The 
temperature  difference  in  the  crude  oil  between  precool  inlet 
and  precool  outlet  xiras  not  very  large  end  was  usually  less  than 
10*^F,  In  the  typical,  case,  the  crude  oil  vas  found  to  increase 
in  temperature  from  two  to  ttiree  degrees  in  travelling  through 
the  pump,  decreased  to  inlet  test  temperature  through  the  precool, 
increased  in  the  test  line  from  a  fraction  of  a  degree  to  a 
maxim.um  of  about  tv/o  degrees.  The  tested  crude  vras  then  blended 
with  residual  crude  of  higher  temperature  in  the  measuring 
tank.  The  solid  lines  plotted  on  Figure  16  show  the  pattern 
obtained  by  using  the  by-pass  after  precooling.  The  lowest 
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test  temperatures  obtained  by  this  method  was  about  30 °F. 

Where  oil  was  by-passed  before  precooling  (Case  II),  similar 
temperature  patterns  were  obtained  from  storage  tank  to  pump  and 
from  test  inlet  to  measuring  tank.  Dissimilar  thermal  beha.vlour 
was  due  to  much  lower  rentes  of  flow  in  the  precool  lines  as 
compared  to  rates  of  flow  in  Case  I.  The  low  rates  of  flow 
caused  slightly  higher  temperatures  in  flow  lines  leading  to 
the  precool  (a  direct  result  of  constant  re-circulation,  without 
cooling,  of  the  bulk  of  the  pump  inlet  stream  to  storage  tank), 
but  in  the  precool  lines,  the  cooling  obtained  v/as  much  greater 
than  in  the  previous  case.  Runs  illustrating  thermal  behaviour 
under  these  conditions,  are  shown  by  the  dash  lines  in  Figure  16, 
Test  temperatures  as  low  as  20°^.  were  obtained  by  this  method. 

B.  Apparent  Viscosity  Characteristics  of  Crude  Oils  Tested 

The  overall  results  of  the  esqDerimental  tests  on  the  crude 
oils  by  the  pilot  pipeline  maj'-  be  summarized  a.s  foUowss 

(a)  Tests  conducted  on  the  pilot  pipeline  on  crude  oils 
were  characterized  by  high  apparent  shear  rates  at 
the  Tjall.  Shear  rates  varied  over  a  wide  range, 

100  to  100,000  sec“^. 

(b)  Apparent  Reynolds  numbers  ranged  from  6l  to  26,000. 

(c)  Apparent  viscosity  values  varied  only  7dth  temperature, 
regardless  of: 

i  wide  variation  in  apparent  shear  rates 
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ii  difference  in  pipe  diameter 
iii  pattern  of  flow,  laminar  or  turbulent 
iv  variation  of  thermal  history 

Thus  apparent  viscosity  values  obtained  by  the  use  of  the 
pilot  pipeline  indicate  that  the  crude  oil  samples  are  Weirtonian 
in  the  shear  rates  range  studied,  or  that  limiting  values  of 
viscosity  at  infinite  shear  for  a  pseudoplastic  material  have 
been  obtained. 

The  above  results  show  that  viscosity- temperature  relation¬ 
ships  for  crude  oils  with  similar  viscosity  characteristics  can  be 
obtained  by  the  pilot  pipeline  from  a  single  run  using  one  test 
pipe  and  one  p8,ttern  of  flow.  Theimal  effects  were  considered 
negligible. 

The  result  of  the  experimental  tests  is  discussed  for  each 
crude  oil  sample  tested. 

(1)  Redwater  Sample  No.  1 

.  The  experimental  data  obtained  for  Redwater  crude  sample 
no.  1  are  listed  in  Table  IV-B.  Tests  1,  2,  3,  4>  S’  snd  6  were 
tests  conducted  on  the  0.249  inch  internal  diameter  line,  while 
tests  ^  ^  v/ere  tests  on  the  0.386  inch  internal  diameter 

line.  Calculated  apparent  viscosity  temperature  relationships 
of  Table  IV-B  are  plotted  in  Figure  17  for  the  0.249  inch  diameter 
line  and  on  Figure  18  for  the  tests  on  the  0.386  inch  diameter  line. 

Tests  on  the  0.249  inch  line  were  run  for  pressure  drop 
ranges  of  12  psi.,  22  psi.,  51  psi.,  58  psi.,  70  psi.,  97  psi., 

180  psi.  and  303  psi.  The  Reynolds  number  variation  in  this 
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range  was  from  12,950  to  96  while  apparent  shear  rates  at  the 
wall  varied  from  746  to  100,000  sec.~^.  Turbulent  flow  data 
could  not  be  obtained  below  47®F,  because  of  inadequate 
refrigeration  capacity.  Laminar  data  could  not  be  obtained 
in  the  higher  teniperature  ranges  (70  to  80°F.)  since  the 
corresponding  pressure  drops  were  too  low. 

The  solid  curve  (Figure  17)  can  be  divided  into  three 
regions.  The  first  region  from  80  to  67°F.,  where  reliable 
values  \7ere  obtained  for  turbulent  flow  conditions,  67  to  47°F, 
where  apparent  viscosity  data  was  obtained  for  both  laminar  and 
turbulent  flo\v  and  in  the  region  47°F,  or  lower,  where  laminar 
flow  conditions  prevailed. 

The  dash  lines  indicaite  apparent  viscosities  computed  for 
conditions  of  transitional  or  uncertain  flow.  Two  sets  of  values 
were  computed  for  transitional  flow,  one  set  assuming  laminar 
flow,  the  other,  turbulent  flow.  Similar  curves  T^ere  obtained 
for  each  set  of  transiti.on  data  and  as  either  laminar  or  turbulent 
flow  became  more  definitely  established,  the  dash  curves  approached 
the  solid  line  defined  by  data  taken  under  definitely  laminar  or 
turbulent  conditions  of  flow. 

The  computed  data  plotted  on  arithmetic  paper  gave  a  curve 
that  was  nearly  linear  with  temperature  in  the  range  80  to  64°F. 
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Below  64.°F.  apparent  viscosity  increased  rapidly  with  decrease 
in  temperature.  The  computed  points,  for  established  patterns, 
also  gave  only  one  curve  and  the  majority  of  data  lay  within 
experimental  error  of  this  curve.  The  single  line  relationship 
found  between  temperature  and  apparent  viscosity  indicate  a 
Newtonian  material,  but  due  to  the  extreme  apparent  shear  rates 
at  the  wall,  pseudoplasticitj?-  cannot  be  discounted.  The  values 
obtained  by  the  pilot  pipeline  may  represent  viscosity  at  infinite 
shear,  hence  the  single  curve  relation  between  temperature  and 
apparent  viscosity. 

The  region  of  transition  from  near  linear  relationship 
between  apparent  viscosity  and  temperature  to  the  steeply  curved 
relationshJ.p  was  very  indefinite.  If  this  region  of  change  of 
curvature  represents  the  beginning  of  was  crystallization,  waxing 
can  be  assumed  to  begin  in  the  range  60  to  65*^F. 

The  computed  data  of  tests  in  the  0.386  inch  line  are  shown 
in  Figure  18.  The  apparent  shear  rate  va.ried  in  these  tests  from 
1470  to  10,300  sec.  Again,  the  data  defined  a  single  curve 
indicating  Newtonian  behaviour  isdthin  the  e:^qperimental  shear  range. 

Comparison  of  tests  on  the  0.386  and  0.249  inch  lines  show 
almost  identical  curves  for  both  pipes.  Near  linear  relationship 
of  temperature-apparent  viscosity  for  the  0.386  inch  pipe  prevails 
to  about  63°F.  In  the  temperature  range  80  to  35°F.,  the  deviation 
of  the  curve  for  the  0.386  inch  line  is  less  than  3  percent  of  the 
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curve  for  the  0.249  inch  line.  This  deviation  increased  at  lower 
temperatures  to  about  6%  at  25°F.  Thus,  r.dtliin  experimental  error, 
values  can  be  assumed  identical  and  pipe  size  disregarded  as  a 
variable  affecting  the  value  of  apparent  viscosity. 

The  temperature  data  taken  during  the  runs  show  a  wide 
variation  in  the  storage  tank  temperatures  of  the  crude  oil  under 
test.  No  effect  of  this  variation  was  observed  in  the  values  of 
apparent  viscosity.  Thus  thermal  histor^^  effects  upon  the  a^pparent 
viscosity  of  the  crude  seem  negligible  for  the  operating  conditions 
of  the  pilot  pipeline. 

(2)  Redwater  Sample  No.  2 

Red^mter  sample  no.  2  had  a  lower  A.P.I,  gravity  than  sample 
no.  1,  34.2  as  compared  to  34»5»  Experimental  data  were  taken 
on  the  0.249  inch  line  and  are  tabiilated  in  Table  V-B  of  Appendix 
B.  Computed  results  are  given  in  Table  III-C  of  Appendix  C,  A 
plot  of  the  daia  is  shown  in  Figure  19.  i'he  data  compared  favorably 
with  those  of  the  first  Redwater  sample  and  from  80  to  45°F.,  the 
apparent  viscosities  were  rd-thin  5  percent  of  the  average  of  the 
first  tests  taken  on  the  0.249  and  0.386  inch  lines,  (refer  to 
Table  VIII-C) .  At  temperatures  below  45*^F,  the  deviation  increases 
to  a  maximum  of  about  20  percent  at  20°F.  Behaviour  of  the  crude 
indicates  wax  precipitation  below  a  temperature  of  approximately 
66°F.  A  comparison  of  the  tvjo  samples  is  shown  in  Figure  24. 
Although  slightly  higher  apparent  viscosity  values  could  be 
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expected  for  sample  no.  2  from  the  specific  gravity  observations, 
this  crude  shows  lower  viscosity  values  in  the  region  30  to  55°F» 
No  attempt  was  made  to  explain  this  phenomenon.  Shear  rates 
varied  from  1,630  to  28,800  sec~^, 

(3)  Stettler  D2-D3  Blend 

Experimental  data  (Table  VI-B)  was  taken  on  the  0.249  inch 
line  under  conditions  of  laminar  flow  over  the  temperature  range 
of  20  to  80*^F,  The  results  (Table  IV-C)  ha.ve  been  plotted  as 
shown  in  Figure  20.  Shear  rates  encountered,  ranged  from  6,640 
to  1,260  sec. 

This  was  the  heaviest  of  the  crudes  tested  (29.1  A.P.I)  and 
as  shovm  on  Figure  24>  gave  h-igher  apparent  viscosity  values  than 
any  other  crude  oil  tested.  Experimental  results  ga,ve  a  well 
defined  curve.  Near  linear  apparent  viscosity- temperature 
relationship  was  found  to  about  59°F.  Below  59°F  the  apparent 
viscosity  of  the  Stettler  D2-D3  blend  increased  rapidiy  with 
decrease  in  temperature. 

(4)  Duhamel-MaJjno-New  Norway  Blend 

The  A.P.I.  gravity  of  this  stock  was  37.1.  Table  VII-B  and 
Table  V-C  show  the  experimental  and  computed  data  respectively. 
Reynolds  numbers  encountered  (242-25,100),  show  that  data  were 
taken  under  botli  laminar  and  turbulent  conditions  of  flow.  Range 
of  shear  rates  was  from  1,0 60  to  158,000  sec.  apparent 
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viscosity-temperature  relationship  obtained  is  shown  in  Figure 

21.  Near  linear  relationship  between  viscosity  and  temperature 

o 

exists  to  about  58  F.  Apparent  viscosity  ranges  from  0.00224 
at  80°f,  to  0.01535  lb.  per  ft.- sec.  at  20°F. 

(5)  Leduc-Woodbend  Blend 

Tables  of  ejq:erimental  and  computed  data  for  '38.5  A.P.I. 

Leduc-Woodbend  stock  are  Table  VIII-B  and  Table  VI-C,  respectively. 

Figure  22  shows  the  plot  of  computed  data.  Waxing  seems  to  hs,ve 
o 

started  near  50  F.  Range  of  shear  rates  during  these  tests  were 
1,080  to  192,000  sec.  Laminar  and  turbulent  conditions  of 
flow  were  used  to  obtain  data.  The  apparent  viscosity- temperature 
curve  for  this  stock  was  similar  to  that  of  Duhamel-Malmo-New 
Norway  stock  except  at  the  lower  temperatures  (Figure  24),  where 
the  effect  of  temperature  on  the  apparent  viscosity  of  the  Leduc- 
Woodbend  blend  was  much  less  than  for  the  Buhamel-Malm.o-New  Norway 
blend. 

(6)  Nisku  Blend 

Experimenta.1  data  for  the  32.6  A.P.I.  Nisku  blend  is  given 
in  Table  IX-B  and  computed  data  in  Table  VII-C. 

Most  of  the  data  were  taken  under  conditions  of  laminar  flow. 
Plotted  data  (Figure  23)  shoivs  that  wax  precipitation  may  have 
started  at  about  58‘^F.  Below  5S°F  apparent  viscosity  increased 
very  rapidly  with  decrease  in  temper.'"  ture.  Range  of  shear  rate 
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encoimtered  in  laminar  flow  tos  825  to  4,950  sec.~^.  The  only- 
turbulent  flow  data  obtained  was  at  a  Reynolds  number  of  7,4J20 
and  a  shear  i^ate  at  the  wall  of  50,800  sec,~^^  Apparent  -viscosity 
values  ranged  from  0.0044  lb.  per  ft.'sec.  at  80°F.  to  0.0430  lb. 
per  ft.- sec.  at  20 

C.  Generalized  Effect  of  Crude  Oil  Gravity  on 
Apparent  Viscosity. 

Figure  27  shows  a  generalization  relating  the  A.P.I.  gravity 
at  60°F.  to  the  apparent  viscosity  and  the  tempera-bure  of  the 
crude  oil.  The  curves  show  a  family  of  lines  and  correlation 
is  good  to  35°F»  -^-t  temperatures  lower  than  35°F.  considerable 

deviation  taJces  -place.  A  plot  of  this  nature  can  be  used  to 
approidjHc  te  the  viscosity-tempera-ture  relationship  in  pipeline 
computatio'ns  for  crude  oils  in  the  absence  of  actual  test  data. 
Such  a  plot,  however,  should  be  used  only  for  similar  base  crudes 
and  crudes  suspected  of  anomalous  behaviour  should  be  tested. 
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CONCLUSION 

The  results  of  tests  on  the  crude  oils  by  means  of  the 
pilot  pipeline  gave  a  single  apparent  viscosity-temperature 
relationship  for  each  crude  oil.  Factors  such  as  pattern  of 
flow,  pipe  diameter,  rate  of  shear  and  thermal  history  were 
found  to  have  negligible  effect  on  the  apparent  viscosity. 

The  crude  oils  tested  indicated  Ne’vrbonaan  behaviour  within 
the  apparent  rate  of  shear  range  studied.  These  observations 
seem  to  contradict  the  findings  of  Sidjak  (l7)  and  Bauer 
and  Huston  (5)^  however,  the  results  of  these  previous 
investigators  were  taken  from  small  laboratory  viscosimeters 
in  which  shear  rate  values  were  very  small.  Maximum  shear 
rates  reported  were  15  reciprocal  seconds  by  SidjaJc  and  40 
reciprocal  seconds  by  Bauer  and  Huston,  Apparent  shear  rates 
at  the  pipe  wall  for  the  pipeline  viscometer  were  usually  of 
the  order  of  thousands  of  reciprocal  seconds.  Hence  possible 
pseudoplastic  nature  of  the  crude  oil  cannot  be  disregarded. 
Since  rate  of  shear  and  pipe  diameter  had  no  effect  on  the 
apparent  viscosity,  values  obtained  could  be  limiting 
viscosity  values  at  infinite  shear  for  a  pseudoplastic 
material . 

The  change  in  slope  of  the  viscosity- temperature  curve 
from  the  near  linear  relationship  found  at  higher  temperatures 
may  indicate  the  beginnings  of  wax  precipitation.  Such  zones 
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of  change  in  slope  of  the  viscosity- temperature  curve  were  in 

o 

the  range  58  to  67  F,  for  most  crudes. 

At  the  request  of  Interprovincial  Pipe  Line  Co.  Ltd.  the 
pilot  pipeline  was  constructed  to  duplicate  the  Reynolds  numbers 
encountered  in  actual  crude  oil  pipelines.  Reynolds  numbers  in 
the  pilot  unit  ranged  from  6l  to  26,000  and  apparent  rs-tes  of 
shear  at  the  pipe  wall  ranged  from  746  to  192,000  reciprocal 
seconds.  In  duplicating  the  high  Reynolds  numbers,  extreme  rates 
of  shear  have  been  introduced  and  for  a  pseudoplastic  material, 
the  majority  or  all  viscositj?"  data  obtained  from  the  pilot 
pipeline  would  be  comparable  to  viscosity  values  at  infinite 
shear. 

Although  turbulent  flow  at  the  higher  Reynolds  numbers  could 
be  obtained,  the  refrigeration  was  found  to  be  inadequate  in 
cooling  the  oil  at  the  higher  flow  rates.  As  a  result,  a  large 
portion  of  the  experimental  data  was  taken  under  laminar  flow 
conditions. 

In  the  discussion  of  refrigeration  studies  in  Appendix  D, 
it  is  shovn  that  the  capacity  of  the  present  refrigeration 
unit  will  not  permit  low  temperature- high  flow  rate  pilot 
pipeline  studies.  The  need  for  further  insulation  at  the 
flanges  is  pointed  out.  Before  further  research  is  continued, 
the  inadequacies  of  the  refrigeration  unit  should  be  corrected. 
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APPENDIX  A 


SAMPLE  COI/IPUTATIONS  OF: 

(1)  EXPEEIMENTAL  DATA 

(2)  FOEMULAE  USED  IN  CO]^IPUTATION  OF  HESULTS 

(3)  DETAILED  RESULTS 

(4)  VISCOSITY  CONVERSIONS 


I  ' 


A  xiaie™  -,. 


’  »y 


'  'k 


V>'f‘ 


10  aHOITATITMDO  ajiM»a  ., 

•.1  ’ 

ATAa  JATMMIfla^Xa  (x)‘ 


8TJV3IH  'io  i.omru'Moo  in  aaau  sajumhot  (&)) 


.A  i:J 


3TJ08SH  (Ci^IATaa  (c) 

i:!oieHa?ioo?rri8oo8iv  (a) 


lA. 


!•  Sample  Computations  of  Experimental  Data  Reported  in  Tables 
IV-B  to  IX-B  (Appendix  B) . 

Average  Test  Temperature 


Test  gauges  provide  for  the  measurement  of  temperatures  at 
points  a  and  e  as  shoijm  by  the  above  diagram,  however  T^^,  the 
average  test  temperature  of  the  test  section  between  points  b 
and  c  is  required. 

Assuming  a  linear  temperature  distribution  along  the  whole 
length  a-e  and  letting  T  be  the  change  in  temperature,  ( T^^^^  -  T^) 
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=  Tfc  -  (AT)  test 

(3-a) 

(AT)  test  =  19-5  (AT) 

(4-a) 

29.63 

T'b  =  T^  .  1.25  (AT) 

(5-a) 

29.63 

from  equations  (2-a),  (3-a.)  and  (4-a)s 
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from  equation  ( 5-a) : 


^av  *"  "^out 


(AT)- 


29.63 


whence  ££  (T„„+-Ti„) 

59.2 


2  X  29.63 
out~^in' 


{AT) 


(6"0t) 


Sample  Computations; 

A  test  in  Stettler  D2-D3  crude  oil  gave  the  following 
experimental  data  (0.249”  line). 

Upstream  pressure^  =30.5  psig. 

Downstream  pressure*"  =2.8  psig. 

Test  inlet  temperature*^  =  72.6  F. 

Test  exit  temperature*"  =  73.0°F, 

Volume  of  crude  collected  =  10.40  gal. 

Time  for  collection  -  8.053  min. 

Temperature  of  oil  in  receiver  =  72°F, 

Temperature  of  oil  in  storage  -  75°^* 
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Temperature  of  oil  at  piimp  exit  -  76°F* 
Temperature  of  F-12  bath  =  65°F. 


^Denotes  average  of  three  readings. 

Average  test  temperature; 

Test  in  Test  out 

Observed  temp.  72.6^F.  73.0°F. 

Corrected  temp.  73.2°F.  73.6^F. 

By  equation  C6-a)  • 


\v  =  73.6  -  (0.4)  =  73.5°F. 

59.6 

Pressure  drop  across  test  section; 
Upstream  pressure  =  30.5  psig. 
Downstream  pressure  -  2.8  psig. 


Pressure  drop  r  27.7  psig. 

Flow  rate; 

Volume  collected  =  10.40  gal. 
Time  of  collection  =;  8.053  min. 


Rate  of  flow  =  10*40  ^  1,291 

-  '  min. 


8.053 

This  amount  viras  collected  in  the  receiver  at  a  temperature  of 
72°F.  The  rate  of  flow  was  corrected  to  the  flow  temperature  of 
73.5°F; 

Density  of  crude  at  72°F.  s  54*59  Ib./ft.^ 

Density  of  crude  at  73. 5°^*  =  54.55  Ib./ft.^ 
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Rate  of  flow  at  test  tempera lure  =  54,59  y-  291 

54.55 

Z  1.292  gal/min. 

or  rate  of  flow  3  1.292  x  -  0.2071  ft^/ 

min  6.24  gal 

2,  Sample  Computations:  Forijulae  Used  in  Compiling  the  Detailed 
Results  R.eported  in  Table  I-C  of  Appendix  C. 

Velocity: 

Flow  velocities  were  computed  by  the  equation: 

V  =  Q  (7-a) 

A 


where  V  =  velocity  in  ft/sec. 

Q  =  flow  rate  in  ft.  /sec, 

2 

A  r  Internal  cross-sectional  aarea  of  pipe,  ft. 

Equation  (7-a)  may  be  simplified  for  velocity  computations  in 

each  pipe.  For  the  pipe  of  0,249  inch  internal  diameter, 
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where  Q*  -  flow  rate  in  ft  /min, 
Friction  factor: 

The  Fanning  equation  is: 
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where  f  =  friction  factor 
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=  pressure  drop,  Ib/ft. 

0  z  density,  Ib/ft. 

V  =  velocity,  ft/sec. 

L  z  length  of  section,  ft. 

D  =  diameter,  ft. 

Since  lengths  and  diameters  are  constant  for  any  pipe,  a 
simplified  equation  for  friction  factor  for  the  0.249  in.  pipe  is! 
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where  (^P) '  =  pressure  drop  in  lb/ in. 

Apparent  Reynolds  Number 

In  the  theory,  it  was  pointed  out  that  friction  factor  could 
be  evaluated  by  the  use  of  apparent  Reynolds  numbers.  In  laminar 
floyj,  the  combination  of  the  Fanning  equation  and  the  modified 
Poiseuille  equation  gave; 
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In  turbulent  flow  a  modification  of  the  Heltzel  equation 
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The  Heltael  equation  wa.s  chosen  for  friction  factor  calculations 
on  the  basis  of  Chmilar’ s  (9)  work  on  the  pilot  pipeline. 

Friction  factors  were  calculated  and  plotted  against  the 
Reynolds  number  on  logarithmic  scales.  The  range  of  Reynolds 
numbers  investigated  was  10,000  to  55>000,  The  Kp^rman-Prandtl 
line  (equation  12)  and  the  Heltzel  line  (equation  14)  were 
transposed  on  the  graph.  By  this  means,  it  was  found  that  the 
Heltzel  equation  best  represented  the  pilot  pipeline  data. 
Apparent  viscosity  (lajninar) : 

Apparent  viscosities  are  computed  by  the  modified  Poiseuille 
equation: 

_  32/^  LV 

o^/a  =  ^  Sc 
^  32  LV 

For  the  o»249  inch  line,  the  equation  reduces  to: 

/a  =  ^  32,2^  _  (3,21  x  10“^)  n 

32  X  19,5  V  V  ft, sec. 

Apparent  viscosity  (turbulent): 

The  apparent  viscosity  in  turbulent  flow  is  derived  from 
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For  the  0.249  inch  pipe,  the  equation  becomes 

/a  -  VP  z  0.0208  VP  lb. 

T\)  ^  ft.  sec. 


(13-a) 


Apparent  Sheer  rate  at  the  v/all: 

The  shear  rate  at  the  vrall  was  shown  by  equation  (20)  to  be; 
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This  equation  is  valid  for  the  laminar  flow  of  a  Mewtonian  fluid 
and  in  turbulent  flow  where  a  laminar  sublayer  exists  next  to 
the  wall.  The  same  equation,  modified  by  the  substitution  of 
apparent  viscosity  is  used  to  compute  the  apparent  shear  rate 
at  the  wall  for  the  laminar  and  turbulent  flow  o|^  non-Newtonian 
fluids.  The  apparent  shear  rate  equations  were  only  used  to 
compare  the  effect  of  a  shear  rate  tem  on  the  apparent  viscosity 
and  are  not  intended  as  expressions  of  true  shear  rate  at  the  wall. 

Substituting  apparent  viscosity  for  viscosity,  the  apparent 
rate  of  shear  for  a  non-Nefidionian  is: 
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In  laminar  flow,  the  modified  Poiseuille  equation  is 
substituted  in  the  apparent  shear  rate  equation; 
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or  since  V  =;  4 
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(15-a) 


For  the  pipe  of  0.249  inch  diameter,  the  equation  reduces  to: 


^  X  103sec-l  (16-a) 

where  Q’  3  flow  rate  in  ft^/min. 

In  turbulent  flow,  the  apparent  shear  rate  at  the  wall 
(l4-a)  for  the  0.249  inch  pipe  is: 


(^)a  r  Q».2.49,  ^  144  .^.32^  (4Z')=  A  ■  (l7-a) 

12  X  4  X  19.5  jUa 

2 

where  A P  is  in  Ib/in. 

3.  Sample  Computations  of  Detailed  Results  Reported  in  Tables  II-C 
to  VII-C  In  cl.  of  -^i-ppendix  C 


(a)  Laminar  Flow 

Using  the  experimental  data  of  Section  1,  Appendix  A: 

Average  test  temperature  =  73*5^I'« 

Pressure  drop  =  27.7  psig. 

Flow  rate  -  1.292  gal ./min. 

Z  0.2071  ftVmin. 

=  54.59  Ib/ft.^ 
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Velocity; 


From  equation  (S-a) , 

V  r  (  Q»  X  10^)  ft 

2.03  sec 

V  =  0.2071  X  10^  =  10.20  ft/sec. 

2.03 


Friction  factor: 

Frem  equation  (9-a), 

f  n  9.8S  (AP)» 
V^|0 


f  -  9.S8  X  27.7 

(10.20)^  C54.59) 


=  4.S25  X  10 


Apparent  Reynolds  number: 
From  equation  (lO-a) , 

(%),  =  i4 
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0.0483 


=  1330 


Apparent  viscosity: 

From  equation  (l2-a) 

/t/a.  =  3.21  X  10"^  (^P«  ) 

'  V 

=  3.21  X  10  ^  ,21*7  .  _  0.00372  lb 
^10.20^  - 

Apparent  shear  rate  at  wall: 

From  equation  (l6-a) , 
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(b)  Turbulent  flow 

A  test  on  Redv;ater  sample  no.  1  gave  the  following  experimental 
data  (0.249  inch  line): 

Average  test  temperature  -  77.5°F. 

p 

Pressure  drop  z  303  Ib/in 

Flow  rate  =  0.S90  ftVmin. 

Density  r  52.3  ft/sec. 

Velocity: 

V  =  ^  X  10^  =  0.390  X  10^  =43.8  ft/sec. 

2.03  2.03 


Friction  factor: 

f  =  9.38  =  9.38  X  303 

(43.8)2  ^2.8 

Apparent  Eeynolds  number: 
by  equation  (3J.-a) 


=  0.0295 


(N  )  =  (Q.v^ii.) 

IX  a  £ 


3.77 


3  77 

=  (0.364\  ’  =  12,950 

0.0295 

Apparent  viscosity: 

by  equation  (13-a) > 

/a  ss  0.0203  Vp  -  0.0208  x  43.8  x  .52.3 

/  -jt)  12,950 

K  a 


=  0.00371  lb 


ft. sec. 


Apparent  shear  rate  at  wall: 
by  equation  (l7*-a). 
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4.  Sample  Computations  for  Viscosity  Conversions  Reported  in 
Table  X-C,  Appendix  C* 

Kinematic  Viscosity: 

The  kinematic  viscosity  of  a  fluid  is  defined  by  the  equation: 

(18-a) 


where  -  viscosity  in  Ib/ft-sec, 

^  =  density  in  lb/ ft.^ 

k.  =  kinematic  viscosity  in  ft^ 


sec. 


From  Table  IX-C,  the  average  apparent  viscosity  of  Redm^ater 

crude  oil,  sample  no.  1  is  19.72  x  10”^  Ib/ft.-see,  at  20°C,,  the 

density  at  this  temperature  is  53.93  Ib/ft. 

Kineniatic  viscosity,  k  -  19.72  x  10~^ 

53.93 


m  3.66  X  10  ^  ft^/sec. 

Saybolt  Universal  viscosity: 

The  relationship  of  Saybolt  Universal  viscosity  to  kinematic 

viscosity  is  given  in  Perry  (15)  as: 

k=  (A9  -  2.  )  cm2  (stokes)  (l9-a) 
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where  A  =  0.0022 
B  =  1.80 

0  =  efflux  time  in  sec  (Saybolt  Universal) 

This  equation  is  valid  for  other  effltix-type  viscometers, 
with  change  in  values  of  the  constants  A  and  B.  Substituting 
apparent  viscosity  in  equation  (l9-a)  one  obtains  - 

1.80 

=y9(0.0022  -  9  )  poise  (20-a) 

converting  to  Eaglish  units: 

1.8 

/*a  -  6.72  p(0.0022  9  -  _  )  lb  (21-a  ) 

'  9  ft-sec 

Equation  (21-a)  is  most  convenient  in  the  fonn 

k  -  6.72  (0.00229  -  1*^  )  ft^  (22-a) 

sec.  ■' 

where  the  values  of  kinematic  viscosity  (k)  can  be  computed  for 
selected  efflux  times  (9) .  Plots  of  this  relationship  are 
available  and  values  of  Saybolt  Universal  seconds  given  in  Table 
X-C  of  Appendix  C  are  taken  from  such  a  plot  for  the  corresponding 
kinematic  viscosity. 
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APPENDIX  B 


EXPERIMENTAL  DATA  OF  PILOT  PIPELINE. 


•  V 


TABLE  1-B 


Test  Line  Dimensions  of  the  Pilot  Pipeline. 
Pipe  size; 


Nominal  outer  diameter,  J- 


Internal  diameter,  in.  0.249 

Cross-sectional 
area,  ft.^ 

Test  length,  ft.  19.5 


0.311 


5 

K 

0.386 


3.38  xlO""^  5.28x10’’4  8.123ao~'^ 


19.5 


19.5 


T/iBLE  n-B. 


Record  of  Crude  Oils  Tested. 


Crude. 

Date 

Received. 

A.P.I. 

Gravity. 

Redwater,  i 
Sample  No.l. 

Apr.  7/52 

34.5 

Redwater, 
Sample  No. 2. 

Feb. 5/53 

34.2 

Stettler  D2-D3 

Blend 

Feb.23/53 

29.1 

Duhamel-Malmo-  Feb. 26/ 53 
New  Norway  Blend. 

37.1 

Leduc-Woodbend 

Blend 

Apr.  20/53 

38.5 

Nisku  Blend. 

Apr. 22/53 

32.6 

Test 

Numbers. 

Remarks. 

1  to  9 

Received  from 

I.P.L.  by 

M.  Chmilar. 

10,11,12. 

Obtained  from 

I.P.L.  from 

Redwater. 
gathering  line. 

13 

Pipeline  sample 
obtained  by 

Canadian  Gulf  Oil  Co. 

14  to  17 

Sample  obtained  by 
Canadian  Gulf  Oil  Co. 

18 

Obtained  from 

I.P.L. 

19 

Obtained  by  I.P.L. 
at  Nisku  and  is 
mainly  Campbell  crude 
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TABLE  III-B 

Density  Deteiroinations  by  Chain  Gravitometer. 


Crude 


Redwater,  sample  no.l. 


Redwater,  sample  no. 2. 


Stettler  D2-D3  blend. 


Temp.  Density. 


gm/ml. 

Ib/ft? 

39.0 

0.8564 

53.47 

54.0 

0.8534 

53.28 

49.0 

0.8539 

53.31 

55.0 

0.8521 

53.20 

60.0 

O.85O8 

53.12 

67.0 

0.8480 

52.94 

46.0 

0.8559 

53.43 

50.0 

0.3538 

53.30 

87.0 

0.8408 

52.49 

72.0 

0.8466 

52.85 

72.5 

0.8469 

52.87 

64.4 

O.S504 

53.09 

63.6 

0.35^0 

53.13 

60. 8 

0.8519 

53.18 

55.0 

0.3550 

53.38 

49.8 

0.8567 

53./^8 

4110 

0.8592 

53.64 

21.0 

0.8682 

54.20 

73.0 

0.3475 

54.59 

60.4 

0.8780 

54.31 

65.3 

0.8765 

54.72 

68.0 

0.8759 

54.68 

37.4 

0.8865 

55.34 

44.8 

0,8853 

55.27 

50.0 

0.8838 

55.17 

53.6 

0.8811 

55.00 

57.2 

0.8802 

54.95 

60.4 

0.8795 

54.90 

44.6 

0.8844 

55.21 

50.0 

0.8832 

55.14 

54.5 

0.8818 

55.05 

57.2 

0.8808 

54.99 
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Crude 


Temp. 


*F.  ;.>i/ral* 

68.0  0.8340 

68.2  0.8343 

39.2  0.8429 

47.3  0.8418 

50.3  0.8408 

53.2  0.8400 

32.0  0.8474 

33.3  0.8466 

35.9  0.8462 

37.8  0.8455 


Leduc-Woodbend  blend.  71.6  0.8271 

63.0  0.8300 

61.2  0.8304 

46.4  0.8340 

36.0  0.8378 

33.3  0,8385 


74.3  0.8553 

61.7  0.8596 

52.7  0.8627 

50.0  0.8640 

30.2  0.8708 


Duharael-Malmo-New  Norway 
blend. 


Density. 

Ib/ftI 

52.07 

52.09 

52.63 

52.56 

52.49 

52.44 

52.91 

52.86 

52.83 

52.79 


51.64 

51;82 
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52,07 
52,31 

52.35 


53.40 

53.67 

53.86 

53.94 

54.37 


Nisku  blend 
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Table  IV-B 

Experimental  Data  of  Pilot  Pipelirje  Tests  for 
Redwater  Sample  No.  1. 


Tests 

in  0.249 

inch  line. 

1  to  6 

inclusive. 

Tests 

in  0.386 

inch  line. 

7  to  9 

inclusive. 

Temperatures  in  deg.F. 

Pressure 

drop 

1. 

storage. 

►  pump 

receiver. 

shell, 

.  test 

Ib/in? 

exit 

average 

76 

76 

76 

76 

76.5 

12,00 

78 

80 

75 

67 

68.8 

11.75 

78 

79 

72 

58 

64.5 

11.95 

76 

77 

70 

49 

59.3 

12.18 

76 

76 

65 

38 

52.4 

12.47 

73 

73 

61 

30 

41.8 

12.95 

75 

75 

60 

26 

32.8 

13.10 

75 

74 

5B 

19 

17.9 

13.75 

75 

75 

57 

18 

19.4 

13.55 

82 

82 

58 

22 

43.1 

12.78 

78 

80 

60 

33 

62.3 

11.90 

74 

75 

62 

48 

66.2 

11.85 

76 

76 

67 

68 

73.2 

11.53 

75 

77 

73 

75 

76.0 

20.2 

78 

80 

74 

63 

70.1 

20.9 

77 

79 

73 

52 

66.7 

19.85 

76 

76 

71 

42 

62.6 

19.33 

72 

73 

67 

32 

55.3 

20.0 

69 

70 

65 

25 

51.0 

20.5 

67 

68 

60 

20 

45.2 

21.0 

65 

66 

59 

18 

40,9 

61 

63 

56 

13 

35.7 

22.0 

59 

59 

54 

11 

31.0 

22.9 

56 

57 

50 

9 

26.7 

23.3 

47 

50 

62 

11 

22.0 

23.8 

56 

60 

54 

23 

22.6 

23.5 

56 

59 

51 

41 

39.5 

22.2 

57 

60 

48 

44 

44.3 

21.3 

57 

59 

52 

30 

51.1 

21.2 

59 

60 

55 

54 

58.8 

20.5 

65 

66 

58 

68 

67.5 

20.4 

66. 

67 

62 

72.5 

21.8 

Flow 

rate. 

ft^/min. 

0.1675 

0.1600 

0.1510 

0.1393 

0.1220 

0.0905 

0.0538 

0.0394 

0.0402 

0.0950 

0.1440 

0.1540 

0.1700 

0.1970 

0.215 

0.224 

0.222 

Q.200 

0.1S73 

0.1667 

0.1530 

0.1366 

0.1200 

0.1075 

0.0942 

0.1073 

0.1523 

0.1680 

0.213 

0.237 

0.252 

0.236 
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Test  Temperal.ures  in  deg.F. 

No* 


Pressure  ELow 
drop.  rate. 


3 


4 

5 


I 

i, 


I 

6 


I 


storage. 

pump 

Receiver. 

Shell . 

test 

Ib/in? 

ft^min 

75 

e^t . 

75 

73 

avern  ge 

75.8 

50.7 

0.318 

73 

75 

72 

30 

62.8 

52,6 

0.314 

69 

70 

67 

21 

57.5 

53.7 

0.317 

63 

64 

64 

18 

51.5 

54.0 

0.320 

56 

57 

65 

15 

45.2 

48.7 

0.351 

51 

53 

61 

11 

41.4 

49.7 

0.353 

56 

58 

61 

12 

46.7 

51.9 

0.336 

59 

60 

59 

15 

54.2 

55.0 

0.318 

62 

64 

58 

22 

61. 6 

53.4 

0.317 

65 

66 

60 

38 

65.2 

52.6 

0.316 

68 

69 

62 

65 

70.9 

51.6 

0.316 

70 

72 

65 

73 

74.5 

51.3 

0.317 

78 

80 

72 

37 

71.5 

57.4 

0.340 

74 

77 

70 

25 

66.6 

58.4 

0.338 

74 

76 

68 

58 

71.5 

68.0 

0.370 

76 

78 

68 

51 

65.5 

68.2 

0.364 

71 

72 

65 

35 

61.0 

69.2 

0.364 

68 

69 

64 

26 

56.4 

70.5 

0.365 

61 

63 

62 

22 

50.2 

73.7 

0.370 

60 

15 

45.8 

74.5 

0.374 

47 

53 

55 

12 

38.6 

67.2 

o.a6 

44 

45 

50 

11 

35.7 

68.8 

0.408 

a 

43 

48 

8 

32.8 

70.5 

0.397 

44 

45 

51 

9 

35.1 

70.1 

0.416 

45 

46 

48 

10 

40.4 

78;  80 

0.388 

48 

50 

47 

19 

45.3 

80.4 

0.383 

51 

52 

48 

29 

50,4 

79.2 

0.383 

55 

56 

50 

44 

55.3 

76.9 

0.378 

60 

61 

54 

58 

61.8 

73.4 

0.378 

63 

64 

58 

64 

65.8 

72.3 

0.375 

66 

68 

60 

69 

69.7 

71.1 

0.380 

72 

80 

65 

58 

77.5 

303 

0.890 

81 

85 
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59 

75.8 

165 

0.624 

77 
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70 

52 
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169 
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72 

75 
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42 

65.4 

174 
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65 
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32 
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20 
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16 
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Test  Temperatures  in  deg.F. 

No. 


Pressure  Flow 
drop  rate. 


Storage. 

pump 

receiver. 

shell 

test 

Ib/in? 

min 

exit. 

average 

61 

62 

58 

6 

49.7 

44.7 

0.301 

60 

62 

56 

6 

53.6 

190 

0.625 

59 

60 

55 

6 

51.6 

70.1 

0.359 

59 

61 

55 

5 

52.6 

101 

0.437 

55 

57 

54 

6 

45.0 

44.3 

0.328 

57 

60 

54 

6 

49.0 

190 

0.621 

56 

57 

53 

6 

46.3 

73.5 

0.370 

55 

57 

51 

7 

48.5 

103 

0.439 

54 

58 

57 

9 

43^4 

48.6 

0.344 

56 

60 

55 

9 

47.8 

192 

0.618 

57 

60 

55 

9 

49.2 

146 

0.535 

77 

78 

73 

65 

78.7 

19.50 

0.605 

77 

78 

74 

65 

75.5 

19.55 

0.602 

73 

75 

74 

51 

69.8 

19.85 

0.599 

71 

72 

71 

42 

64.7 

20.0 

0.598 

67 

67 

66 

34 

59.8 

20.45 

0.594 

63 

62 

64 

30 

55.8 

20.8 

0.^3 

57 

58 

63 

17 

50.7 

22.1 

0.604 

5S 

60 

67 

13 

56.3 

21.0 

0.594 

62 

65 

65 

18 

60.3 

20.4 

0.591 

67 

68 

65 

34 

65.7 

19.93 

0.591 

71 

73 

66 

71.2 

19.72 

0.597 

76 

78 

71 

75 

75.9 

19.57 

0.604 

76 

75 

76 

61 

73.0 

10.28 

0.420 

71 

72 

72 

59 

66.7 
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65 

66 

65 

47 

61.0 
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60 
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38 

55.3 

9.26 
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58 
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55 
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Test 

Temper 

atures  in 

deg.F. 

Pressure 

Flow 

No. 

drop 

rate. 

storage. 

pump 

receiver. 

shell 

test 

Ib/int 

ft'Vmin’ 

exit. 

average 

9 

70 

71 

70 

69 

71.0 

18.82 

0.582 

68 

69 

69 

63 
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TiiBLE  VI-B 


Experimental  Data  of  Pilot  Pipeline  Tests 
for  Stettler  D2  -  D3  BlCUiid  in  0.249  Inch  Line. 
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TABLE  VIII-B, 


Experimental  Data  of  Pilot  Pipeline  Test 

.  for 

Leduc-Woodbend 

Blend  in  0 

.249  Inch  Line. 

Test 

Temperature 

in  deg.F, 

Pressure 

Flow 

No. 

drop 

rate 

storage  pump 

receiver 

shell 

test 

/  2 
lb. /in. 

ft.  Vmin 

exit 

average 

18 

71 

78 

72 

59 

72.7 

323.9 

1.020 
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76 

72 

67 

73.8 

3.^68 

0,03663 
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4.016 
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466.9 
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TABLE  IX-B. 

Experimental 

Data  of  Pilot  Pipeline 

Test 

for  Nisku  Blend  in  0, 

,249  Inch  Line. 

Test 

Temperature 

1  in  deg. 

F. 

Pressure 

Flow 

Wo. 

drop 

rate 

storage  pump 

receiver 

shell  test 

lb,/in.^ 

ft.^/min 

exit 

average 

19 

72 

76 

74 

58  73.6 

207.4 

0.6803 
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72 

68 
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18.88 

0.2214 
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19.34 

0.2092 
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Test  Temperatures  in  deg. 
No. 
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pump 
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48 
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Flow 
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rate 
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ft. '/min 
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APPMDIX  C 


DETAILED  RESULTS  OF  COJilPUTATIONS  OF  PILOT 


PIPELINE  TESTS 


1C 


TABLE  I-C. 


Summary  of  foniiulae  Used  in  Computations 


Reported  in  Section 

2  (Appendix  A) 

• 

Internal  pipe  diameter,  in. 

0.249 

0.311 

0.386 

Velocity,  ft/sec. 

Q’xlO^ 

Q'xlO^ 

Q'xLO^ 

2.03 

3.17 

4.87 

Friction  factor. 

9.88AP' 

''V 

12.35AP> 

15.32^P‘ 

V'^/o 

Apparent  Reynolds  number. 

laminar 

f 

64 

f 

M 

f 

turbulent 

f 

Apparent  viscosity,  Ib/ft-sec 

laminar 

-3 

3.21x10(ap«) 

(  V  ) 

-3 

4.99x10  (XpO 

T'v) 

7.71x10 

(  V  ) 

turbulent 

0,0208(VP} 

0.0259(Vp) 

(Wt,) 

s. 

0.0322  (Vp) 

(“ba 

Apparent  shear  rate  at  wall. 

~1 

sec. 

5. 10(^0^ 

laminar. 

I8.92fc0^ 

9.740^0^ 

turbulent 

1.231  (AP*) 
(/a) 

1.541  (AP') 

(/b 

1.912  (AP') 
(/a) 

Note;  Units  of  syrnbols  used  in  the  above  equations  are;- 

Q'  =  flow  rate,  ft^/min. 

AP*  =  pressure  drop,  Ibs/in. 

=  density,  Ib/ft.^ 

V  -  velocity,  ft/sec. 


aaold'JBj’i/qftloO  xtx  baaU  oeXo&nol  "io  ’TiB  roiuje 


.  (A  xLtbnsqqA)  S  lioitooQ  as  fasdnoqefl 


b&t.O 

ne.o 

«iix  t'xad’araflib  ©q±q  iBffxei’nl 

Sj&c'P 

^OIx’P 

^Oix'P 

.  0  sa  Vi  ole  V 

'  vi.e 

co.s 

<\^v 

‘^88,9 

A 

.loc^ojal  nojtJ^oli^ 

,*i©cf£DLUfx  abXoirceH  ixitnsqqA 

1 

iBolml 

1 

1 

1 

{  1  ^ 

&aBSjj(Sisj& 

.osaAlVl  ^Tjd-jtaooaJrv  d-na^aqqA 

^jOJicIV.V 
(  V  ) 

(  vT 

e- 

( »qA)0IxXS.t 
(  V  ) 

xanXmaX 

(qy)  ssco.o 

(qV)9eSQB0 

Wwi 

3 

(qV)80S0.0 

ineXuQ'xud’ 

^0L^0I,5 

( 

^"".oas 

,Ilsw  d’B  &&ert  xaexla  ,tii8XBqqA 

^oij^p^v.e 

^OlxJpS^.BX 

.XBnlxneX 

SI9,I 

(b\) 

'A’ 

OtM  xes.i 

w 

d-fielwJxud 

-iOTtfl  cjnolJ-£j;pa  avodB  al  bsaif  alodat^a  lo  ailnU  zecfoM 

.nlinX^iM  tioII  =  ’P 

^/.fljtXerfl  ,qoib  euxreaortq  =  ’‘iA 
^.d-lVX  »^ilanab  =  ^ 

.oaaV^  ^^d’loolav  r  V 

i 


Detailed  Results  of  Computations  for  Pilot  Pipeline  Tests 
on  Redwater  Sample  No,  1* 
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Detailed  Results  of  Computations  of  Pilot  Pipeline 
Tests  on  Duhamel-Malmo  -  New  Norway  Blend. 
(Uncertain  flow  results  omitted  from  graph) 
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Detailed  Results  of  Cociputations  for  Pilot  Pipeline 
Tests  on  Nisku  Blend. 

(Uncertain  flo?;  results  omitted  from  graph) 
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Apparent  Viscosity  Comparisons  for  Redwater 
Crudes.  (Values  from  Figures  17,  18  and  19) 
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APPENDIX  D 


REFRIGERATION  STUDIES 


ID 
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Refrigeration  Studies  Conducted  on  the  Pilot  Pipeline 


To  determine  whether  more  insulation,  larger  surface  area  of 
heat  transfer  in  the  precool  bath,  an  auxiliary  precool  refrigeration 
unit  or  a  refrigeration  unit  of  greater  capacity  is  required  to  give 
the  desired  lovirering  in  oil  temperature  at  the  higher  rates  of  flow,  an 
investigation  was  made  of  the  refrigeration  capacity  of  the  existing 
vmit  at  various  rates  of  crude  oil  flow  through  the  0.386  inch  internal 
diameter  test  line.  The  flow  rates  used  were; 

Test  I  Zero  flow  through  test  line,  maximum  flow  through 
by-pass  (18,36  g.p.m.) 

Test  II  Small  flow  through  test  line  (5.64  g.p.m.)  bulk  of 
flow  through  by-pass. 

Test  III  Medium  flow  through  test  line  (8.35  g.p.m.)  moderate 
flow  through  by-pass. 

Test  IV  Large  flow  through  test  line  (12.71  g.p.m.),  small 
flow  tirrough  by-pass. 


For  each  case,  the  by-pass  valve  was  adjusted  for  the  desired 
rate  of  crude  oil  flow  through  the  test  line,  the  remaining  portion  of 
the  crude  oil  pumped,  was  by-passed  to  storage  without  cooling.  The 
oil  was  all.owed  to  circulate  through  the  storage  tank,  refrigerated 
shell  and  receiver  and  data  of  Table  I-D,  Appendix  D,  taken  after 
equilibrium  conditions  were  obtained. 


lo  aerte  ©OB'i'Tirs  »Tora  'ledieriw  ©aliined-eb  oT 

nolcfB^oa-t'xls'x  loooeiq  rta  ,xld-iid  loooanrq  ed&  rut  islsnsi;)-  ;^BorI 

svi-j  CKt  beTiwpei  8l  \:iJ:oBqxio  lod-Baaa  lo  ^xm  rroliBnsaXilo’i  3  lo  Xlm; 
ns  ^woll  'io  boXot  'larlslrf  orf;}'  d'B  o’mti.'TCdqfaed-  Ilo  id  ^nircswol  bsixaab  eiU 
aid.t-aix9  9iIX  "io  AjXXoBqjao  fiolv+BTOgxTC'^©*!  dilX  lo  ebBffi  Qon  noi&a^iiBbval 
iBrfred’id  do  id  dBt*0  arid  lisirOTdd'  won  Ixo  obimo  lo  sei'&'i  bx/oXiby  Xb  &lnLr 

t  V  ¥ 

:©i9w  besif  aadjB*!  woH  9j1T  ,snll  d-Bod  ledoioiib 
ii§mndd  wdn  irtinnixiaiTi  ^exdI  daeX  dsircndd  woJl  oocoS 

(.m-.q,3  d8,8I)  esaq-Y*^ 

lo  :ILyrf  (.m.q»g  ^idl  d'aad  dguondd-  woll  IlanS 

»3  3Bq"--’id  rigtronrld  woXl 
ed^Yabom  (.ic.q.g  ^C.S)  enxX  daed  dsxr>ii[d  woH ■“ cK/xboIiS 

V  *1 

ffe.r,roarid  woll; 

JCIame  ^(.lE.f'.g  IV«SI)  ordi  deod  tlj>iuo*idd  woH  egiaj 

.  'Bjsq-rcf  d^uorcrid  woll 

bsYisfib  odd  tol  E«>d0jyt,bJ3  aaw  ovJQjv  rianq-^cf  odd  ,00^0  doae  lo*? 

^  f 

"lo  ncxdioq  aidrtxxeio't  odd  tonxl  deed  add  riaiiprrdd  \mll  Uto  obxrio  T;6  adBi 

) 

ariT  ,3111X000  di/oridiw  ogaioda  od  bseenq-Yd  a'Bw  *bsqttti-q  Ilo  ebijnco  add 
bedirtoai’i'loTT  ^irwd  aginoda  add  riamjndd  edBlooilo  od  bawoIXn  airfl  Xlo 
'todlo  nE*yi!!d  ,cr  xXbitaqqA  ^0.-1  eXdjeT  lo  Ajub  bus  tovxsosrr  baa  IXeda 

.beotadcfo  oiaw  saoXdxbnoo  nitiXTxfilXirfje 


I  daeT 

II  daoT 

III  daeT 


VI  dear 


2D 


TABLE  1-D 


Experimental  Data  for  Refrigeration  Studies  on  the  Pilot  Pipeline 

(Refer  to  Figure  1-D) 


Test*  Symbol 

Units 

I 

II 

III 

IV 

Compressor: 

F-12,  suction  press. 

psig. 

1.5 

1.5 

3.0 

3.0 

discharge  press.  P2 

psig. 

95 

115 

120 

122 

discharge  temp.  T2 

°F. 

97 

162 

172 

192 

Power  input 

Kw. 

2.05 

2.08 

2.26 

2.22 

Condenser: 

Water,  inlet  temp.  T^ 

®F. 

60 

60 

62 

62 

exit  temp .  T, 

°f! 

68.6 

76 

84.6 

84.6 

floiY  rate  Ib/sec. 

0.343 

0.351 

0.372 

0.382 

F-12,  exit  temp.  T^ 

Op 

64 

66 

76 

72 

Shell: 

Oil,  precool  inlet  temp.  T, 

Op 

- 

37 

50 

71.5 

test  line  inlet  temp.  T^ 

Op 

- 

29.8 

41.6 

68.0 

test  line  exit  temp.  Tg 

Op 

— 

30.4 

42.7 

69.0 

test  line  inlet  press. Py 

psig. 

- 

39.0 

91.0 

174.0 

test  line  exit  press.  Pg 

psig. 

- 

4.2 

6.9 

15.8 

test  flow  rate  Q.ji 

g.p.m 

0 

5.64 

8.35 

12.71 

by-pass  flow  rate  Qg 

g.p.m 

.  18.86 

- 

- 

- 

Bath  temp.  Tg 

Op. 

1.0 

. 

H 

H 

28.5 

56.5 

Average  vap-^ur  barrier 

surface  temp.  T* 

Up. 

57 

60.1 

60.9 

65.6 

Average  flange  Surface 

temp.  Tp 

"F. 

25 

28 

40 

59 

Room  temp.  Tg 

°F. 

77.1 

78.1 

73.0 

72.5 

Barometric  press. 

psia. 

13.7 

13.7 

13.7 

13.6 

T  0.386  inch  internal  diameter 
w.  Redwater  crude  oil,  34.5  A.P, 

line  used  in  all  tests 
,1.  in  test. 
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Detailed  Outline  of  Computations  Used  in  BefrLgeration  Studies 


A.  By  Heat  Balance  Method, 
(a)  Condensers 


P-12.  L*c^o*d 
64“  F. 


The  amount  of  F-12  circulated 
around  the  condenser; 


where 


Data  (Test  I); 

F-12  enters  condenser  as  vapor, 
temperature,  T2  =  97oF. 
pressure,  -  95  psig. 

F-12  leaves  condenser  as  liquid, 
temperature, , T^  =  64®F* 

Cooling  water  leaves  condenser, 
temperature,  T^  =  68.6®F» 
Cooling  water  enters  condenser, 
temperature,  T^  = 
rate,  =  0.343  Ib./sec. 

Atmospheric  pressure  =  13,7  psia. 


is  found  by  a  simple  heat  balance 

Qc  -  Qw  Cp  (aT)t5^  =  Qf  (  I-d) 

=  T4  -  t^  =  change  in  temperature  of 
cooling  water 

(^H)p  =  enthalpy  change  of  F-12  or  the  enthalpy 
of  F-12  vapor  at  temperature  T2  and 
pressure  P2,  minus  the  saturated  liquid 
entha.lpy  of  F-12  at  temperature  T^. 


J 


ft-  . 

3ejtbxid-3  nl  beuU  'lOOid-Bcf-irqmoO  lo  aulIiJ'jjO  bolixxfea 


:<(I  cJ-ftef)  jeJ'aQ 

'  ■  '  ■  <  ^  At.’’ ^ 

'8B‘ '18eJ3ipbfiCK>  'rVj©j-f£9  Si.~'U 
,:i  »_^T  t^yi£x:'-j8*rsqii'0d' 

.3x30^  ,  ^s'rxiQB'Jiq’ 

tbiirpxJl/SB  ’■toQfSEibado  b^vb^X  SX-i 

'  '  ' 

.‘5°+V^  -  ^^ysrtud'ii'soqi^cf- 
tietjfiejgno^  3ev,  '.sl  T0d'«\'r  ■gnJcIooO 

i  '■ 

Ao^ .  8d  ~  te'a^j8T:0qiflBJ- 

tiaansMoo  ko&sifr  ,3niloo0 

^©lijd-BiegaiaJ' 

oe3\.dX  Zi^K.O  =  ted-Bi 

\  ' 

* -'.ISc  V.^X-=  OTiu^aerfq  oXieriqaonvhA 


.boilddll  soxioIbS  d'aoH  xS  ‘A 
cnsanobnoO  (a) 


% 

It 


<i 

r 

♦ 

- 

'  *  t 

'  -sw 


SIX-'?  lo  cfiitfoaiB  edT 

\ 

:*io3nabflaj  odj  ’basicriB 


aan«  LDd  ,  d^aori  olmnia  b  \(S  b'auoJ  al.  b&&Bli/0’ilb 


(b-I  )  qsP  tP  -  ol> 

I. 

lo  o-Tij^aiaomevt  nX  P^niido  =  -T  -  erterfw 

'  ,t  ■•• 

•SO&  MV  jullooo 

t 

arid  10  ^I-^l  lo  Oafwii>  ‘ViqXflry-n©  = 
brw  laJ  rrwqiBBi  f.i  logjsv  SJi-VlQ 
blu  ii  boXAiiilBii  odd-  Bi/fii  'i  ,;;1  sncxfoaoTcq 
.;iT  0T:wdr.-i’'‘  if>.t  ,t«  IJI-*?  lo'''cqXBdd'no 


> 


5B 


Gp  =  mass  rate  of  flow  of  F-12 
Cp  =  specific  heat  of  water 
and  =  heat  loss  to  the  surroundings  by  the 
condenser. 

The  heat  gained  by  the  cooling  water  (q^)  is: 

Ow  z  Op;  Cp  =  (0.343  (1.0  )  (68.6  -  60)‘>F 

=  2.95  Btu/sec. 

The  enthalpy  (hg)  of  F-12  superheated  vapor  at  97°f  and  108.7 
psia.  from  Perry  (15)  is  87.99  Btu/lb.  The  enthalpy  (hf)  of  F-12  liquid 
at  64^F,  assuming  saturation,  is  22,50  Btu/lb. 

Hence  (aH)^  =  hg  -  hf  =  (87.99  -  22.50)  Btu/lb. 

=  65.49  Btu/lb. 

Assuming  adiabatic  operation,  Pe  =  0,  thus  from  equation  (1-d) : 
2.95  =  65.49  Qp 
or  Qp  =  0.0450  Ib./sec. 

which  is  the  rate  of  F-12  flow  through  the  refrigeration  unit  for  Ce.se  I. 

The  mass  rate  of  flox?  of  F-12  may  also  be  computed  by  considering 
the  volume  displacement  of  the  compressor.  Two  cylinders  of  inch 
diameter  with  a  3#  inch  stroke  are  driven  at  the  rate  of  65O  cycles  per 
minute. 

The  theoretical  volume  displaced  per  min  =  ^  (2-d) 

■  2  ‘ 

where  D  =  diameter  of  cylinder 
S  =  length  of  stroke 


ij 


lo  wort  lo  ©Ntai  yaara  ~  <^3- 

,1  ' 

toi:ya  'to  ifsaorl  oilloeqa  «  qO 
0xii  '^cf  aanioxiuoi'itrti  edf  q&  yaol  taarl  =  bna 

.Toa/T©Lnoo 

:8l  (irP>  'loSsvt  ^ntlooo  erf.^  '^cf  boniueg  i'ssxl  sdT 


.003 Wa  29.2  = 


V.80I  bns  J-B  ffbnav  h&^a&iijequQ  SI-'l  lo  (^d)  xgXcrid-na  ©ilT 
biupij:  C;i-'5  to:  (-jd)  ■’cqXsjU'ro  oriT  .d’lV-sl'a  99.Vfe  al  (:?!)  itmt  .eleq 


.dlV^a  (y?.?.S  3X  ^nol.+eijud-B8  anxuitraaa  d-a 


.cllWa  (O^.SISi  -  :..?8)  »  =  ^(Ha)  ©oneH 

.-fi'wa  9^.2^  = 


:(ci*-X)  noid’Bup©  aixU  ^0  =  qP  ^noi^Biaqo  oid-':cfBJ!:bB  gnlnujoeA 


.oo3\.cfI  OeAO.O 

.T  n'*.aO  To't  t  xnx/  noldl'iErraal'i'ioT  ©dd’  djjpt^d  \foIt  SJC-'l  lo  ©dai  exld  el  rioxxlw 


ax-rebLufoo  bifrJ’fTt'flioo  ed  oeXa  ^bki  'ti-'i  to  voTi  'lo  ed’BX  eajjin  sriT 

'51  '  '  1  ' 

ibat  to  TxelinXIico  ov/X  .1' •  {aoiqnBoo  edd'  Ip  d’xtomooB  Ccji;:?)  oijujXov  ed.t 


.pdixniBt 


’T  ^^jTT  ®  nim  ’T  .  r ''ijqouXqait*  ©cttxlov  Xaoxd'saoarid’  »tlT 
_ 


iftbr.M';;©  to  'ted’oiuMlb  =  Q  ©TioriTr 
.'O'vtdo  lo  <';rtoX.  ~  3 


6D 


Substituting  values  in  equation  v2-d) ; 


2 

theoretical  displacement  =  ^  (2,5)  (3«5)  (650) 

2  (1728) 

=  12.84  ftVmin. 


mm. 


Assuming  F-12  to  enter  the  compressor  as  saturated  vapor, 
(2.490  ft  /lb.),  at  the  suction  pressure  of  1.5  psig.,  or  15.2 
psia.,  the  capacity  of  the  refrigeration  unit  is: 


(2.490  ft Vlb)  (60  sec/min) 

The  flo?7  of  0.0860  Ib/sec.  is  much  higher  than  the  flow  rate 
of  0.0450  Ib/sec.  calculated  by  the  heat  balance  method,  however, 
the  flow  rate  of  O.O86O  Ib/sec.  does  not  account  for  the  volumetric 
efficiency  of  the  compressor  nor  possible  superheating  of  inlet* 
vapors. 


(b)  Shell 


64*  F 

0  045  Ib/sec. 
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Data:  F-12  enters  expansion  line  at  temperature. 


*^5  Z  64°F 


F-12  expands  in  line  at  suction  pressure. 
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=  1.5  psig 

No  oil  flow  through  test  line 
Shell  temperature,  Tg  =  i.q 
Assumptions: 

1.  Saturated  F-12  liquid  enters  expansion  line  at 

2.  Saturated  F-12  vapor  enters  compressor  at  15.2  psia, 

3.  Flow  of  F-12  is  distributed  equally  in  each  line. 

From  Perry  (15) : 

Saturated  liquid  enthalpy  of  F-12  at  64  ^  z  22.50  Btu/lb. 
Saturated  vapor  enthalpy  of  F-12  at  15.2  psia  -  75,84  Btu/lb. 
Therefore  change  in  enthalpy,  (AH)'p,  -  53.34  Btu/lb, 

By  a  simple  heat  baiance  around  the  shell,  the  heat  gained  by 
F-12  in  evaporating  8.t  15.2  psia  equals  the  heat  transferred 
through  the  shell  plus  the  heat  transfer  from  the  flow  of  crude  oil 
in  the  test  line  under  equilibrium  conditions: 

qg  t  So  =  Qj.  (3.^) 

where  q^  r  heat  transfer  through  shell  plus  losses  in 
refrigeration  capacity. 

=  net  heat  transfer  from  oil,  defined  by  the 
thermodynamic  flow  equation. 

(aH)f  =  change  in  F-12  enthalpy. 

For  Case  I,  where  there  is  no  oil  floT^,  substitution  of  values 
in  equation  (3--d)  gives: 

=  0.0450  (53.34)  BWsec. 

-  2.4  Btu/sec.  or  8,640  Btu/W. 
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Thus  8,640  Btu/lir  is  transferred  into  the  shell  in  maintaining 
the  shell  at  a  temperature  of  1.0°F,  Actually  some  of  this  heat  is 
transferred  through  the  external  piping  and  may  be  considered  as 
losses  in  the  refrigeration  capacity.  No  estimation  of  losses  was 
attempted. 

(c)  Sample  computations  for  case  of  oil  flow  through  test  line. 


F-ie 

0.0739  'Vsec. 
66*  F. 


15.2.  psiQ.. 

J 


29.8  F. 


Shell 


Oil 

5-64  g.p.m. 
37*  F. 


The  heat  transferred  (q^)  from  the  oil  may  be  calculated  by 
the  thermodynamic  flow  equation: 

[aU  +  ^(Pi)  +  t  asg]  -  .  n,  -  q 

2  gq  go-*  ■  ^  ® 

Where:  -  ma^ss  rate  of  flow  of  oil 

-  Jc^dT  -  change  in  internal  energy,  U. 

P  r  Pressure 

0  r  specific  volume 

V  -  velocity 

^  n  height  above  datum 

%  n  net  heat  transferred,  positive 

when  heat  is  transferred  from  oil  stream. 
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z  shaft  work 
“  specific  heat  of  oil 
T  r  temperature 

g  -  acceleration  due  to  gravity 
gf.  z  conversion  factor 

Since  there  is  no  shaft  work  done  on  the  oil  in  the  precool 
and  test  lines,  -  o.  The  change  in  potential  energy  is  small 
and  is  assumed  negligible.  The  change  in  specific  volume  and 
specific  heat  of  the  oil  is  assumed  negligible.  The  equation 
becomes; 

lo  3  -Qj  [o^'ST  +  I  ]  (4.<a) 

2  gc 

Data;  (Test  II) 

o_ 

Oil  enters  precool  at  temperature,  T^^  -  37.0  F 
Oil  enters  test  line  at  pressure,  =  39.0  psig. 

Oil  enters  test  line  at  temperature,  Tr^  -  29,8°F 
Oil  exits  test  line  at  temperature,  Tg  -  30.4°F 
Oil  exits  test  line  at  pressure,  Pg  -  4,2  psig. 

Ba.te  of  flow  of  oil,  Qtji  -  5»'34  gwi /min .. 

Density  of  oil; 

at  37.0°F  -  53.6  Ibs/ft^ 

30.4°F  =  53.7  Ibs/ft^ 

29.8°F  =  53.7  Ibs/ft^ 
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Apparent  viscosity  of  oil; 

at  37.0°F  -  0.0102  Ib/ft.sec. 

30.4°F  =  0.0121  Ib/ft.sec. 

29.B°F  -  0.0123  Ib/ft.sec. 

F-12  enters  e^cpansion  line  at  temperature, 

T5  =  66°F 

F-12  expands  in  line  at  suction  pressure, 

1  -  1.5  psig. 

Calculated  rate  of  flow  of  F-12,  Qp  z  0.0739  Ib/sec. 
Cross-sectional  area  of  test  line  0.000812  ft.'~ 

From  Perry  (15) : 

Specific  heat  of  oil,  C^  _  0.41  Btu/lb.®F. 

Internal  diameter  of  precool  line  -  0.622  inches 

2 

Cross  sectional  area  of  precool  line  -  0.00211  ft. 
Assumptions; 

1.  F-12  enters  e^cpansion  valve  as  saturated  liquid. 

2.  F-12  enters  compressor  as  saturated  vapor. 

From  a  heat  balance  ar  und  the  shell,  the  heat  gained  in 
evaporating  F-12  (qp)>  is  equal  to  sum  of  heat  transferred 
from  the  oil  (q^),  the  heat  transferred  through  the  shell 
plus  the  heat  transferred  through  the  external  piping  (q^),  or 

but  qj.  =  Qp  (AH)J, 

%  =  -  Qt  [c^  (*T)  +  ^(*P)  4-  1 
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Substituting  values  in  equation  (4-d) : 


-  =  (5..-64  X  53.7  x  60  lb)  (0.4I  ^  )  (37.0  -  30.4)  °F 

6.24  hr.  Ib^F 


Z  7880  Btu/hr. 

The  pressure  drop  in  ^(ap),  of  equation  (4-d)  is  the  press>ure 

drop  from  precool  inlet  to  test  outlet.  This  includes  55  feet 

of  ^  inch  precool  pipe,  8  feet  of  0.386  inch  internal  diameter 

calming  section  and  19.5  feet  of  0.386  inch  internal  diameter 

test  section.  Obstructions  to  flow  in  the  precool  section  are 

five  90°  elbows,  three  tees  and  one  sudden  contraction.  Equivalent 

lengths  for  these  obstructions  are  from  Brown  (7): 
o 


5-90  elbows  -  5  x  1.5  ft.  = 

3  -  tees  =  3  X  3.5  ft.  z 

1  sudden  -  0.6  ft.  - 

contraction 


7.5  feet 
10.5  feet 

0,6  feet 


Total  equivalent  length  of  obstiuctions  =  18.6  feet 

Leng-th  of  ^  inch  precool  pipe  -  55*0  feet 

Total  equivalent  length  of  -g  inch  pipe  =  73.6  feet 


The  velocity  of  the  oil  in  the  I  inch  precool  line  is: 


|/  r/5»64  ft.^  \  ^  1. _ 

0,24  X  60  sec.  0,00211 


z  7.14  ft/sec. 


The  ma^dLmum  Reynolds  number  in  the  f  inch  pipe  will  be; 


=  (0.622)  (7.ai4)  (53.7)  =  1950 
/a  ^  (12)  (0.0102) 
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T 
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The  pres'jure  drop  in  laminar  flow  for  the  inch  pipe 
is  calcul.ated  using  the  modified  Poiseuille  equation,  that  is; 


=  32LV^  ^ 

Substituting  values: 


4P  -  ,g2_(73.^6iL7JA), (0.011^)  -  15.3  it/i^ 

(32.2)  (0.622)*’ 


The  pressure  drop  in  the  test  section  is: 

^7  "  ^8  =  ^9.0  -  4.2  -  34.8  Ib/in.^ 

The  pressure  drop  in  the  8  foot  calming  section  is: 

S  X  34.^  lbs,  ~  14.3  Ib/in,^ 

19.5  in^ 


The  total  change  in  pressure  in  the  precool  and  test  lines 
is  -  64.4  psig.  The  product  of  the  specific  volume  and  the 


pressure  drop  becomes: 


“  Qrp'^(^P)  =  f5.64  X  6Q  X  53.7  lb,  ) 

6.24  hr.  S3. 7  lb. 


(64.4 


.4.  X  144  lb. 

ft.' 


>)( 


Btu 

778  lb. ft) 


i:  646  Btu/hr. 

2 

The  kinetic  energy  term,  ^  V  may  be  computed: 

2  gc 

velocity  in  inch  precool  line  3  7.14  ft/sec. 


velocity  in  test  line  -  5.64 _  _ 

8.12  X  10~^-  X  6.24  X  60 

z  18.6  ft. /sec, 

-  ftp  (AL  )  =  _  2920  ,•18.6^  -  7.P  1  _  17.2  Btu/hr. 

2  Sc  9  X  32.2  ’ 


eqlq  do  of  {  -  uld  •xo't  woH  rj^nlmL  ni  qmb  our-  ae-iq  edT 
;ai  divi.t  ,noii'Bwpo  sIIxjjBalo^  ^nJttu  bb^^alMol^o  ai 


VJSif.  - 


.3 


;a8uXBv  SaiiiAtidecfijS 


.nixdx  e.2X  r  i;:-'-.';q-,P).^4,v).(^,.gf),  SS;  :  ^ 
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^.fdtVX  8.Ae  r  -  0.9C  = 


:ax  rroidooa  ani.TiXiiO  doo'i  8  ofld"  ni  qoib  enx/aa&'ic  exIT 


^.nlNpfl  Z.U  r  x  _8 
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ori^  DOF-  e  .“nlov  oilxooqa  erfi  lo  diox/bonq  oriT  .gleq  -  ai 

:B0iii009c[  qonb  OTfiwee'iq 

■  ^ 

■<!s-  7.  :?  •■¥ 


^  .'viX  It  X.id)  (,,  -''-•-2-S-i^^^^)  =  - 


"■'  .nxiWS  - 

c 
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and  the  total  heat  transferred  from  the  oil  is: 

qo  =  (7,880  I  646  -  17.2)  ^  z  8,509  Btu/hr. 

hr 

From  Perry  (l5)! 

hg  at  15.2  psia  =  75.04  Btu./lb. 

at  66°F  -  23.00  Btu/lh. 

hence  -  52.84  Btu./lb. 

Thus  qp  z  0.0739  (52.8)  Btu/sec.  -  3,90  Btu/sec.  or  16,250  BtiVhr, 
Finally  the  heat  transfer  through  the  shell  from  the  surroundings 
is: 

Ig  =  (16,250  -  8509)  3tu/hr. 

-  7,741  Btu/lir. 

In  a  similar  manner,  the  heat  transfer  through  the  shell  was 
computed  for  each  test  and  tabulated  in  Table  II-D,  Appendix  D, 

(d)  Total  Energy  Balance 

From  the  results  of  the  preceding  computations,  an  overall 
energy  balance  may  be  written  around  the  refrigeration  unit. 

Consider  the  system  to  include  the  refrigeration  lines,  compressor, 
motor,  condenser  and  shell.  Then  the  energ;^/’  balance  is  given  by: 

%  =  9s  ^  ®m 

where  q^  is  the  heat  energy  transferred  from  the 
system  by  the  water 

q^  is  the  heat  energy  transferred  into  the 
system  through  the  shell  and  the  external 
piping,  equation  (3-d). 


;al  XI  o  sdd'  pionl  iijema'ia/iifTd'  (tjeerf  iBchot  9ild-  bna 


.’iti\fi&a.  z  (s-vi  -  I  oaStV)  =  qP- 

•xrf 


•  (^l)  ^vio^  morc^' 
•dlVij^Q  ^8,5V  =  siQq  S.5X  iB 
.diVc^a  oo.es  =  ’iPdd  &B  ^ 
.dXX.jjda  ^a*S^  -  ^(H^)  sorted 

,*iri\/uda  0?S,ai  *10  ^osaXudG  r  .oeaVicj’S  (S.S^)  PCVO.O  -  .^p  ajjrlT 

aanJtbntrontrB  ©dd  tncrtl  Xlaria  ed&  dgjLrortdd  islarLsid  dsari  odd  icIJufinJtl 


:ai 


.iri'vud-a  (poes  -  oes^dx)  =  qP' 

.irlWa  .XAVtV  r  . 

'' '  I  }-i'* 

aav  lloda  edd  riguortdd  Ts'iaofxt  daeri  ©rid  tTsimam  xaHjitla  a  nl 

.d  xXbnaqqA  ^d-II  ©IdaT  n±.  badslndsd  boa  dasd  rioa©  lol  boduqnioo 

,  , 

'  ^  i 

»  ■  I 

eonalaa  ^i9ivi  IsdoT  (b) 

i-Xeievo  na  ,afioidj3daqraoo  snlbeoeiq  ©rid  lo  Bdliraarc  erid  ffiorr'i 

•  *•  •  ' 

.dinxr  noxds'ta^lTlo’i  ©rid- brsaorta  neddiiv?  ©d  ^Brn  ©onelad  \a'i©f(e 
^noBasTcqnioo  tEonlf"  noidineali'i©*!  arid  ebjoXonJt  od  medaAca  ©rid  leblanoO 
'  :'^d  ri©v±3  -ii  eonaXad  ’^giane  arid  naitT  .XIaria  bna  taansbnoo  <iodora 

t  O?  +  aP  :  wP 

arid  .rcrj^  beTi^elariB^t  ^jgneria  dsari  add  al  ^  eiailw 
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‘i&daw  erid  •^jd  aiada-'Ca' 
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and 


is  the  heat  energy  transferred  into  the 

system  by  the  oil,  equation  (5-d) 

is  the  electrical  energy  supplied  to  the  system. 


Expaos\on 

Valve 


From  condenser  cooling  water  requirements  (part  a) : 

q^  m  (2.95  Btu/sec.)  (3600  sec,/hr.'>=  10  ,  600  Bfu/hr. 
From  total  energy  balance; 


qg  =  8640  Btu/hr. 

%  =  o 

®in  z  2.05  (0.75)Kw  (3413  Btu/Kw.hr.)  =  5250  Btu/hr. 
(assuming  75^  efficiency  of  motor  and  compressor 
unit) 

%  -  qg  T  qQ  I  -  13,390  BtVhr. 

For  each  test,  the  energy  output  at  the  condenser  may  be 


computed  by  the  total  energ;^^  balance  and  compared  to  the 
sensible  heat  gained  by  the  water  at  the  condenser. 


odnl  berrxQlfias^d-  '^lone  iaeil  odd-  al 
(^>“v)  nojfcd-oxjp©  ^Ii:o  odd  ■\^d  tn&dQy^ 


,trji&p.xB  odd-  oi  boxIaqiLiB  'VS'ieno  IflOjhcdooio  ©xld-  si  bxoi 


.*i-iWa  00^^01  sC.iA.oea  OOd?)  (.oesWa  e^-S)  -  D 

<  “  w* 

;80fielccf  ^^lorre  Latod  <Borrt 

.rcfiV/d-a  CU^8  =  gP 


.xdV«dH  o?s?  =  cue)  wxc^ev.o)  eo.s  z 

■:30iqiiiQO  brio  lodom  lo  anlmimao) 

(d-jtnir 

.triWa  oee,px  -  ^,o  ^  1- 

.d  ’;8u;  lofcittidxiijo  sad  d.-?  di;qdi/o  edd  ,daeu+  rioij©,  -xoH, 

'»dd  od  beicAqiDoo  baa  00 nr, Lad  Ipdod  odd  wi  b&dan'-oo 


‘  ( 
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In  the  compression  of  F-12  vapors,  the  compressor  was  assumed 
to  operate  in  a  reversible  adiabatic  manner  (constant  entropy  process). 

In  actual  practice,  consideration  of  entropies  at  compressor  inlet 
and  exit  conditions  showed  that  the  process  could  not  be  adiabatic. 

In  Case  I,  the  entropies  of  F-12  vapor  for  the  inlet  condition  of 
15,2  psia.  and  -20°F  and  the  exit  condition  of  10S.7  psia  and  97°F, 
were  0.1727  Btu/lb.°F.  and  0.1694  Btu/lb.'^F  respectively.  The  heat 
rejected,  neglecting  irreversibilities  is  approximately; 

^57.  44Q  (-  0.0033)  =-1.65  Btu/ib. 

At  the  F-12  flow  rate  of  0.0450  Ib/sec.,  the  heat  rejected  is: 

0.0450  (3600)  (1.65)  r  267  Btu/hr. 

The  heat  e^ipelled  at  the  condenser  bj"  the  total  energy  baiance 
is  decreased  by  the  amount  of  heat  rejected  by  the  compressor,  that 
is,  the  total  energy  balance  becoraes: 

^-u  =  9s  +  I 

where  =  heat  rejected  at  compressor. 

For  Case  I,  the  heat  expelled  at  the  condenser  becomes  13,620  Btu/hr. 
In  general  heat  effects  at  the  compressor  could  become  quite  large  when 
operating  at  higher  pressures  and  flow  rates. 
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B,  Heat  Transfer  Through  Shell  by  Method  of  Overall  Resistance 
Coefficient 


Data:  as  showi  on  following  sketch; 


Ou  Ci?~ =  3-  2S 

S'tju~/'e7cjz  7e/>9-/7.  =  5^?  V 
CO/PQUGm'E0  VAPO/^  3AR^/£fi 

O  ■  O  /in  .  f/y/cA 

/<  *  //3  eAy.fii:/Ar. 

ff^/PGLAS  /A/SULPT/ON 
X  =  002 Z 


STEEL  CAS//VG 

8-£23>' /h  .  O-D. 

7  L2S />>  Z-O. 
2G.2 


Casing  length  between  flanges: 

Calii'jing  section 

Test  section 

Total  Casing  24  ft.  3f  in. 

o 


6  ft.  5  in. 

17  ft.  lot  in. 


Inner  wall  temperature  estimated  at  2.0  F. 

TheKEo conductivity  values  from  McAdams  (13). 

The  conduction  of  heat  through  several  bodies  in  series  is 
given  in  Perry  (15)  as  follows: 

q=  (6.a) 

Itfp 


aioruftc^aieefl  IliPisvO  lo  borlJ-eM  ^cf  XleriS  dsj:;cmlT  islaaa^  d-BsH  ,a 

dr^iolllecO 

; vvJ  ^jiJ^oIXol  no  mofle  as  zadaQ 


raegnsll  .rieewd-sd  rid^nel  gnlBBO 
.nx  2  .d"!  ^  nold-oea  grttralBO 

.at  oOX  .ctl  VI  '  fioJtd'oea  d'asT 


.nl  JC  *d‘’i  AS  gnianO  Ijad-oT 

o 

.'i  O.S  is  bed'BiiLttae  enxjd'jsneqnjsd’  Ilnw  •xanni 
.  (Cl)  a/nabAoM  oKn'l  eoniBv  ''jd'tvliox/bnooofiTisffT 
al  aehtoa  nl  aeiJbod  Ixrxevee  rianorudd’  daeri  lo  noidoxrfjnoo  etfT 


':BwoIIo’t  BB  (Cl)  vno'^  ni  nevts 
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where  q  -  heat  transferred  through  the  composite  wall 
At  z:  the  sum  of  the  individual  temperature  drops 
across  each  material 

Itj,  =  the  sum  of  the  individual  resistances  ^2  etc. 

The  resistance  R  is  defined  as: 

R  =  1  (7-d) 

KA 

where  1  -  thickness  of  heat  conducting  ma.terial 
K  -  theimal  conductivity  of  material 

=  cross-sectional  area  perpendicular  to  the 
direction  of  heat  flow. 

Substituting  values  and. computing  the  individual  resistances 
per  foot  of  pipe: 


(a)  Steel  casing 


R,  =  l, 


S  S 


(0,50/12)  ft, _ 

(26.2  ^  2~^)  (2,16  ) 


hr, ft.'  “F 
Z  0.000733  hr."F.Ft/Btu. 
where  -  2_Qg  jjjgan  area  per  foot  of  pipe 


ft. 


=  -  Dj;)  _  'jr(l.o)  ftVft. 


In  D, 


% 


12  In  8.625 
7.625 


_  2,16  ft.  /ft, 
(b)  fiberglas  insulation 


(S-d) 


r  ( — - )  hr.  ft. 

KjAj  ^  0.022  X  2.71  Btu. 

r  2.53  hr.  °F.  ft./Btu. 
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By  equa,tion  (S-d)  the  log  mean 

area  _  IT  (l2,25  -  8.625) 

12  In  12.25 
8.625 

=  2.71  ft.Vft. 


where  outer  diameter  of  insulation  -  3,25  (l2)  in. 


'fr 


=  12,25  in. 


end  inner  diameter  of  insulation  -  8.625  in. 
(c)  Aluminum  vapor  barrier 


~  =/0.0Q45 _ 

KA  ^12  X  118  X  3.25 

CL  CL 


-  0.000,00098  hr.  ft. 

Btu, 


■)  hr,  ^F.  ft. 
Btu. 


whence  =  (0.000733  ^  2.53  )  hr.  °F.  ft. 

Btu. 

=  2.53  -hr.  ^F.ft. 

Btu, 


q  =  (57.0  -  2.0)‘^f 

2.53  hr.^^F^.  (24.29  ft) 

Btu. 

z  527  Btu/hr. 


The  heat  transfer  through  the  end  flanges  were  computed 
on  a  log  mean  area  based  on  the  area  exposed  to  the  cooling  bath 
and  the  outside  area  of  the  flange  face: 

^  ^  ~ 
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where  -  outer  :.rea 
r  inner  area 


Letting  q  be  the  heat  transfer. 


hence  A^  =  -  7.63^) 

TTx  la  In 


q  z  ka  ^t  - 

X  ~ 


(26.2  X  Q.6OI  X  23)  Btu 


2.25/12  hr. 


Ijn  - 


7.63' 

=  0.601  ft.^ 


z  1930  ^ 

hr.  '  -  - - ^ 

The  unit  is  provided  with  two  exposed  dead  end  flanges 
and  two  partially  insulated  inner  flanges  (Fig.  8),  Sstimating 
the  heat  transfer  through  the  insulated  flange  to  be  one-third 
of  the  exposed  flange,  the  total  heat  transfer  through  the 
shell'^estimated  as  5?  677  Btu/ hr. 

Heat  transfer  through  casing  527  Btu/hr. 

Heat  transfer  through  exposed 


flanges. 


2  X  1930  = 


3860  Btu/hr. 


Heat  transfer  through  insulated 


flanges. 


2  X  1930  = 
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Total 


1290  Btu/hr. 
5677  Btu/hr. 
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The  heat  transfer  by  the  resistance  coefficient  method 
was  computed  for  each  case  and  is  presented  in  Table  II-D, 


TABLE  II-D 

Results  of  Heat  Transfer  Computations. 
Test  Conditions 


Test 

Oil  flow 

Shell 

Precool 

Test 

Test 

Comore 

ssor 

Rate 

Temp. 

Inlet 

Inlet 

Outlet 

Suction 

Temp. 

Temp. 

Temp, 

Pressure. 

gal/min. 

0 

F 

°F 

^f 

Psia. 

Psig. 

I 

0 

1.0 

— 

15.2 

1.5 

II 

5.64 

11.5 

37.0 

29.8 

30,4 

15.2 

1.5 

III 

8.35 

24.0 

50.0 

41.6 

42.7 

16.7 

3.0 

iv 

12.71 

56.0 

71.5 

68.0 

69.0 

16,6 

3.0 

F-12  flow  Rate 

Total  Refrigeration  Capacity 

(1)  By  Heat 

(2)  By  Compressor 

(3)  Estimated 

( 4)  Rated  © 

Balance 

Displacement 

590  R.P.M. 

lb./ sec. 

lb,/ sec. 

Btu/hr, 

9 

Btu/hr. 

1 

0.0450 

O.OS6O 

8,640 

33,400 

©21,0  psig. 

II 

0.0739 

0.0860 

16,250 

27,700 

©15.2  psig 

III 

0.115 

0.0940 

21,400 

23,200 

©10.7  psig 

IV 

0.107 

0.0934 

20,400 

Heat 

Transfer  from  Oil, 

Btu/hr. 

(9)  Heat 

Transfer  from 

Shell  by  Heat  Balance 


(5)  Internal  (6)  Heat  of  (?)  Kinetic  (8)  Total  Method,  Btia/hr, 


^ergy 

Change 

Friction 

Energy 

Change 

(including  losses) 

I  0 

0 

0 

0 

8,640 

II  7,880 

646 

-19 

8,510 

7,740 

III  13,000 

2,360 

-55 

15,310 

6,090 

IV  6, 640 

6,540 

-194 

12,990 

7,ao 
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Heat  Transfer  by  Resistance  Heat  Expelled  at  Condenser, 

Coefficient  Method,  Btu/  hr.  Btu/lir. 


(10) 

Casing  (ll) 

Flanges 

(12)  Total 

(13)  By  Temp,  of  (14) 
Cooling  Water 

By  Total 

Energy 

Balance 

I 

527 

5,150 

5677 

10,600 

13,890 

II 

457 

3,470 

3927 

20,200 

21,570 

III 

302 

2,330 

2632 

30,200 

27,170 

IV 

87 

794 

881 

31,100 

26,070 

Discussion  of  Results  of  Refrigeration  Studies 

Columns  (8)  ajid  (l2)  of  Table  II-D  give  the  estimated  heat 
transfer  through  the  shell  by  heat  ba.3.ance  method  and  by  resistance 
coefficient  method,  respectively.  Comparison  of  values  obtained, 
show  that  estimations  by  the  latter  method  are  lower.  Heat 
transfer  values  sho^wi  in  column  (8)  should  be  decreased  by 
the  amount  of  heat  transfer  due  to  external  piping  and  the 
amount  of  heat  transfer  at  the  refrigerant  expansion  valves. 

No  a.ttempt  was  made  to  evaluate  the  heat  transfer  from  these 
sources. 

In  general  the  heat  transfer  rates  through  the  shell  are  high 
and  show  the  need  of  further  insulation.  Transfer  rates  tlirough  the 
flanges  (column  11),  account  for  sJjnost  all  the  total  heat  transfer 
through  the  shell  and  indicate  the  necessity''  for  adequate  insulation 
of  these  flanges. 

Estimated  refrigeration  capacities  (column  3)  indics.te  the  need 
for  a  larger  refrigeration  -unit,  for  example,  if  the  minimum  storage 
taaperature  is  estimated  at  30°F  and  the  crude  oil  is  to  be  tested  at 
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20°F  at  a  high  flow  rate,  say  12.71  g.p.m.  as  in  Test  IV,  the 
sensible  heat  transfer  from  the  oil  alone  is: 


-,12.71  ft.^ 
6.24  min. 


.1-)  (10°F.) 


=  27,700  Btu/hr. 

The  figures  for  test  IV,  column  (6)  show  that  approximately  65^0 
Btu/'hr  of  heat  will  be  generated  due  to  friction  at  this  flov;  rate. 

The  total  amount  of  heat  transfer  Tuthout  considering  losses  through 
the  external  piping  system  exceeds  the  refrigeration  capacity  of  the 
unit.  An  alternative,  to  the  replacement  of  the  present  unit  may 
be  through  installation  =  of  an  aujdJ.iary  precooling  system.  The 
cooling  capacity  of  such  a  unit  must  be  at  least  equal  to  the 
difference  between  the  sum  of  heat  transfer  through  the  shell 
and  heat  transfer  fi*om  the  oil,  less  the  capacity  of  the  present 
unit. 

Comparing  the  estimated  refidgeration  capacity  (column  3)  with 
the  rated  capacity  (column  4)  >  the  refrigeration  unit  is  operating 
near  rated  capacity  for  the  suction  pressure  used.  More  efficient 
use  could  be  made  of  the  available  refrigeration  if  a  closer 
approach  of  test  temperature  to  bath  temperature  could  be  realized. 

This  would  necessitate  increased  area  of  precool  piping.  Heat 
e>rpelled  by  the  condenser,  colimms  13  and  14j  check  T.ithin 
expectations  and  shomr  that  the  results  of  computations  by  heat 
balance  method  are  fairly  accurate. 
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The  flow  rates  computed  by  the  sensible  heat  gained  by  the 
water  for  Case  111  and  IV  are  considerably  lower  than  those 
calculated  by  the  compressor  displacement  method.  The  only 
explanation  for  this  discrepancy  is  that  inlet  water  temperatures 
measured  by  the  ’pjrrocon’  gauge  may  have  been  in  error. 


